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ABSTRACT 


Development  of  an  effective  method  for  evaluation  of  alternative  energy  sources  in 
the  automotive  industry  has  always  been  a  necessity  for  cost  efficient  design  analysis.  One 
viable  alternative  energy  source  is  electricity.  In  the  present  day  environment  of  shrinking 
fossil  fuel  supplies  and  environmental  awareness,  electric  powered  vehicles  are  becoming  a 
low  cost,  non-polluting,  alternative  means  of  transportation.  The  analysis  of  reliable 
electric  propulsion  can  be  expensive  without  a  modeling  tool  for  evaluating  design 
strategies  before  vehicle  construction. 

This  thesis  explores  electricity  as  an  alternative  energy  source  for  the  automobile 
of  tomorrow.  Under  the  guidelines  of  the  World  Solar  Challenge,  a  solar  powered  electric 
vehicle,  using  a  permanent-magnet  brushless  dc  motor  has  been  modeled  and  simulated  in 
Simulink  (Dynamic  System  Simulation  Software).  The  simulations  were  performed  with 
the  goal  of  determining  the  optimum  configuration  to  efficiently  utilize  the  power  supplied 
from  the  solar  array,  batteries,  and  motor.  The  simulated  vehicle  was  "driven"  over 
various  terrain's  and  at  various  speeds.  The  results  obtained  confirm  this  simulation  as  an 
efficient  design  tool  and  present  an  example  of  an  optimum  vehicle  speed  for  one 
particular  vehicle  configuration. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Development  of  an  effective  method  for  evaluation  of  alternative  energy  sources  in 
the  automotive  industry  has  always  been  a  necessity  for  cost  efficient  design  analysis.  One 
viable  alternative  energy  source  is  electricity.  In  the  present  day  environment  of  shrinking 
fossil  fuel  supplies  and  environmental  awareness,  electric  powered  vehicles  are  becoming  a 
low  cost,  non-polluting,  alternative  means  of  transportation.  The  analysis  of  reliable 
electric  propulsion  can  be  expensive  without  a  modeling  tool  for  evaluating  design 
strategies  before  vehicle  construction. 

The  World  Solar  Challenge  is  a  race  that  was  designed  to  demonstrate  current 
technology  in  the  field  of  Solar  Powered  Electric  Vehicles  (SPEV).  Since  it's  inception  in 
1987,  it  has  provided  industrial  and  educational  "teams"  an  opportunity  to  present  the 
latest  technological  trends  and  breakthroughs  in  solar  powered  electric  vehicle  research  to 
the  public.  The  primary  means  of  solar  powered  vehicle  development  has  been 
characterized  by  theoretical  evaluation  of  the  vehicles  velocity  based  on  four  variables: 
aerodynamic  drag,  rolling  resistance,  mass,  and  power  output  Most  of  the  vehicles 
entered  in  the  World  Solar  Challenge  are  evaluated  with  these  four  variables  before 
construction. 

A  mote  effective  design  tool  would  take  these  variables  and  incorporate  them  into 
a  dynamic  vehicle  model.  Coupled  with  the  motor,  battery,  and  solar  array,  the  design 
engineer  could  then  base  the  final  vehicle  design  on  data  accumulated  from  simulations  of 
this  dynamic  model.  A  Solar  Powered  Electric  Vehicle  Simulator  is  therefore  absolutely 
necessary  as  a  primary  design  tool  for  cost  effective  electric  vehicle  research  and 
development 

B.  THESIS  OVERVIEW 

The  primary  means  for  development  of  solar  powered  vehicle  designs  has  been  to 
devise  a  mathematical  model  based  on  the  performance  variables  of  rolling  resistance,  air 
drag,  gravity,  velocity  and  required  system  energy.  This  model  is  used  to  calculate  the 
power  required  to  maintain  a  continuous  speed  under  steady  state  conditions.  TTiis 
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approach  is  adequate  for  design  evaluation  when  not  considering  transitions  in  terrain, 
speed,  or  cloud  cover. 

The  goal  of  this  study  is  to  develop  a  dynamic  simulation  tool  for  evaluation  of  a 
solar  powered  electric  vehicle  model  based  on  the  system  variables.  These  system 
variables  will  include  the  vehicle's  mechanical  and  electrical  components.  This  approach 
will  allow  the  designer  to  specify  and  change  the  vehicle's  characteristics  in  order  to 
optimize  the  design.  In  final  form,  this  modeling  tool  will  realistically  simulate  a  solar 
powered  electric  vehicle  with  changing  terrain,  speed,  and  solar  array  power  while 
providing  continuous  on-line  data  display. 

Figure  1-1  shows  a  block  diagram  of  the  vehicle  design  considered  in  this  study. 
The  primary  components  are  solar  array,  battery,  voltage/current  bus,  motor,  drive  train, 
and  vehicle.  Models  for  the  electrical  and  mechanical  system  elements  are  presented 
independently  and  tied  together  to  produce  the  system  model. 

1.  Background  of  the  World  Solar  Challenge  and  Solar  Powered  Vehicles 

This  section  provides  a  brief  history  of  solar  powered  electric  vehicles  and  the 

events  surrounding  the  first  World  Solar  Challenge  [Ref.  1,  pp.  1-6].  Discussed  are  the 
basic  rules  governing  vehicle  construction  and  some  of  the  vehicle  designs  of  previous 
races. 

2.  Modeling  the  Vehicles  Electrical  Components 

Discussed  here  are  the  solar  array,  battery,  voltage/current  control,  and  motor. 
The  array  and  battery  characteristics  are  developed  from  experimental  data  and  placed  into 
matrices  which  are  used  as  look-up  tables  in  the  system  model.  The  current  bus  ties 
together  the  array  current,  battery  current,  and  the  regenerative  current  from  the  motor. 
The  output  current  form  the  bus  is  integrated  to  monitor  the  battery  state  of  charge 
(SOC).  The  vehicle  system  model  utilizes  a  permanent-magnet  brushless  dc  motor  which 
is  presented  in  state  space  form. 

3.  Modeling  the  Vehicles  Mechanical  Components 

The  system  variables  which  describe  the  drive  train  and  torque  loading 
characteristics  of  the  vehicle  are  developed  in  this  section.  The  forces  on  the  vehicle 
consideted  to  be  of  primary  importance  are  rolling  resistance,  air  drag,  acceleration  force, 
and  MU  gravitational  force.  The  coupling  equations  wMch  describe  the  drive  train  are  also 
presented  in  this  section. 
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VEHICLE  DYNAMICS 


4.  Simulation  Software  and  Model  Development 

Ease  of  use,  speed,  and  output  capabilities  are  factors  considered  important  when 
choosing  the  simulation  software.  The  primary  software  tools  were  Matlab  and  SimuHnk 
from  MathWorks  [Ref.  2].  Both  provide  for  excellent  modeling  of  the  linear  and  non¬ 
linear  differential  equations  which  form  the  basic  structure  of  the  system.  Each  vehicle 
component  is  developed  and  presented  as  a  Simulink  submodel.  These  submodels  are  then 
tied  together  to  form  the  dynamic  system  of  the  solar  powered  electric  vehicle. 

5.  Model  Response  and  Data  Analysis 

The  system  model  is  simulated  and  analyzed  for  varied  driving  conditions.  Plots  of 
the  responses  are  presented  and  discussed.  The  simulation  results  are  compared  to  the 
theoretical  expected  outputs.  Then,  the  simulation  is  used  to  present  a  basic  driving 
strategy  to  optimize  the  chosen  vehicle  configuration. 

6.  Conclusions  and  Future  Work 

Finally,  conclusions  about  the  usefulness  of  the  solar  powered  electric  vehicle 
simulator  are  presented.  Suggestions  for  expanding  the  model  to  include  simulation  of 
accessories,  power  losses,  thermal  characteristics,  and  external  data  accumulation 
capabilities  are  presented. 
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n.  BACKGROUND  OF  THE  WORLD  SOLAR  CHALLENGE  AND  SOLAR 

POWERED  VEHICLES 


A.  THE  WORLD  SOLAR  CHALLENGE 

1.  World  Solar  Challenge  Race  Background 

The  "World  Solar  Challenge"  race  was  the  idea  and  creation  of  Hans  Tholstrup, 
the  first  man  to  biuld  and  drive  a  solar  powered  electric  vehicle  2500  miles  across 
Australia  [Ref.  1,  p.  1-6].  He  embarked  upon  that  adventure  primarily  to  call  attention  to 
the  need  for  a  reliable  alternative  energy  source.  The  World  Solar  Challenge  was  an 
opportunity  for  private,  corporate,  and  educational  teams  to  prove  their  various  vehicle 
designs  and  give  the  world  an  opportunity  to  view  the  technological  advances  being  made 
in  alternative  energy  sources  for  vehicles  of  the  near  future. 

The  first  World  Solar  Challenge  brought  a  field  of  24  entrants  from  around  the 
world.  It  began  on  November  1,  1987  and  was  won  by  the  General  Motors  team  entry 
Sunracer.  The  average  speed  of  the  GM  Sunracer  was  41.8  mph,  which  allowed  the  team 
to  cover  the  1868  mile  distance  from  Darwin  to  Adelaide  in  44.9  hours. 

The  design  strategies  used  in  the  1987  race  set  the  precedent  for  the  next  World 
Solar  Challenge  race  held  in  1990,  which  brought  36  entrants  to  Australia.  The  course 
was  again  the  1868  mile  Stuart  Highway  between  Darwin  and  Adelaide.  The  race  was 
won  by  the  Biel  Engineering  School  team  entry  Spirit  of  Biel  II.  The  Spirit  of  Biel  II 
averaged  40.5  mph  and  completed  the  race  46.1  hours.  Table  2-1  shows  the  vast 
improvement  in  the  average  results  achieved  by  all  entrants  who  completed  the  race  in 
1987  and  1990. 

2.  Rules  Governing  Vehicle  Design 

The  general  rules  governing  solar  powered  race  vehicle  configuration  can  be  used 
as  a  basic  outline  for  development  of  a  modeling  tool  or  simulation  program.  The  vehicle 
limits  considered  important  in  this  study  are: 

•  maximum  solar  collector  array  dimensions  are  2  m  wide  by  4  m  long  by  1.6  m 
high  -  this  limits  the  array  power 

•  maximum  vehicle  dimensions  cannot  be  greater  than  6  m  long  by  2  m  wide  1 .6 
m  high  +  1%  allowance 

•  minimum  height  of  the  car  is  1  m 

•  minimum  weight  of  driver  is  80  kg 
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These  limits  provide  a  basis  for  design  of  the  submodels  developed  in  the  following 
chapters.  Since  solar  powered  electric  vehicles  are  historically  designed  to  be  lightweight, 
the  need  to  optimize  the  power  delivered  from  the  solar  array  is  paramount. 


Average  Speed,  mph 

Average  Finishing  Time, 
h 

1987 

19 

124.4 

1990 

28 

69.8 

Table  2-1.  Comparison  of  World  Solar  Challenge  Race  averages. 
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m.  MODELING  THE  VEHICLES  ELECTRICAL  COMPONENTS 


A.  THE  SOLAR  ARRAY  COMPONENT 

The  data  for  a  solar  array  output  to  an  ideal  constant  voltage  bus  was  taken  from 
data  collection  experiments  performed  at  Rose-Hulman  Institute  of  Technology  [Ref.  3]. 
The  data  collection  was  completed  by  testing  three  types  of  solar  array  configurations;  flat 
array,  curved  array,  and  flat  tracking  array.  The  following  assumptions  and  conditions 
were  present  during  data  collection: 

•  the  size  of  the  array  conforms  to  the  World  Solar  Challenge  race  rules  which 
state  the  array  must  fit  into  a  box  1.6  by  2  by  4  meters. 

•  racing  is  done  between  9  am  and  6  pm. 

•  the  day  of  solar  collection  is  clear  and  sunny. 

Open  circuit  voltages  and  short  circuit  currents  were  collected  every  30  minutes, 
and  the  available  power  was  then  calculated  for  each  array  configiuation.  Since  the 
simulation  uses  the  current  output  from  an  array  for  a  specific  time  of  day,  the  data  may  be 
programmed  into  a  look-up  table  for  use  in  the  simulation  routine.  Figure  3-1  presents  the 
normalized  current  output  from  the  three  array  configurations  versus  the  time  of  day. 
Although  this  data  is  not  exact  it  does  present  an  approximation  that  is  quite  close  to  the 
expected  results  as  compared  to  the  array  output  of  the  vehicles  entered  in  the  1990 
World  Solar  Challenge  [Ref.  1,  App.  4].  Figure  3-2  shows  the  block  diagram 
implementation  of  the  current  supplied  from  the  solar  array  to  the  vehicle  system. 

B.  THE  BATTERY  AND  BUS  COMPONENTS 

The  battery  component  of  the  simulation  model  is  presented  in  the  same  form  as 
the  solar  array.  Battery  discharge  data  supplied  by  the  manufacturer  [Ref.  1,  pp.  73-75]  is 
placed  into  a  matrix  and  then  incorporated  into  a  look-up  table  for  use  in  the  simulation. 
The  bus  combines  the  solar  array  and  battery  for  output  to  the  motor,  and  takes 
regenerative  current  from  the  motor  to  the  battery  for  charging. 

1.  Modeling  the  Batteries  and  Bus 

Two  types  of  batteries  are  used  in  the  model.  Lead  Acid  and  Silver  Zinc.  Figure  3- 
3  shows  the  battery  voltage  versus  percent  charge  or  state  of  charge  (SOC).  The  output 
voltage  to  the  motor  will  have  a  maximum  value  of  100  volts,  although  the  input  to  the 
voltage  bus  may  be  larger. 
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Figure  3-1.  Available  Current  from  Solar  Array  vs.  Time  of  Day. 
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Figure  3-2.  Block  diagram  of  current  supplied  from  the  solar  array  to  the  vehicle 

system. 
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The  battery  State-of-Charge  model  (SOC(t))  is  developed  from  the  Martin  Model 
[Ref.  4,  p.  87].  It  is  calculated  as  a  function  of  time  from  the  current  drawn  by  the  motor, 


SuaeofOutrge(%) 

-^LeadAcid  -+-SUverSnc 


Figure  3-3.  Battery  Voltage  vs.  State  of  Charge. 


current  supplied  by  the  array,  and  current  regenerated  by  the  motor.  An  expression  for 
this  relationship  is 


SOC{t)  =  SOC(O)  -  100 


lo(  ^mot  ha  ) 

Qo 


(3.1) 


where  is  defined  as  the  battery’s  maximum  current  capacity  in  amp-hours,  4^  is  the 
current  supplied  from  the  solar  array  as  defined  in  the  previous  section,  and  equals  the 
motor  current.  The  motor  current  is  positive  when  the  motor  is  drawing  current  and 
negative  when  the  motor  is  generating  current.  This  regenerative  action  of  the  motor  is 
caused  by  slowing,  stopping,  or  going  down  a  hill.  SOC(O)  is  the  battery's  initial  state  of 
charge  (in  percent). 
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The  block  diagram  shown  in  Figure  3-4  presents  the  implementation  of  the  battery  and  the 
voltage  bus  to  be  used  in  the  overall  vehicle  model. 


Figure  3*4.  Block  diagram  of  battery  and  bus  model  to  be  simulated. 


C.  THE  PERMANENT-MAGNET  BRUSHLESS  DHtECT  CURRENT  (PMBDC) 
MACHINE 

A  simplified  approach  to  development  of  a  permanent-magnet  brushless  dc 
(PMBDC)  motor  wiU  be  presented  from  the  Krause  method  [Ref.  5,  pp.  273-328]  of 
machine  analysis.  This  established  theory  of  PMBDC  machines  wiU  present  the  machine 
dynamic  characteristics  in  terms  of  linear  differential  equations,  which  wiU  then  be  given  in 
state  space  form.  The  matrix  of  state  equations  wUl  be  the  basis  for  construction  of  the 
time-domain  block  diagram  for  the  dynamic  model. 

1.  The  Machine  Volts^e  and  Torque  Equations 

The  two-pole  two-phase  permanent-magnet  dc  machine,  as  shown  in  Figure  3-5, 
wiU  be  used  to  develop  the  machine  equations.  In  the  circuit  diagram,  represents  the 
identically  distributed  as  and  bs  windings  on  the  stator  and  the  winding  resistance.  The 
d  axis  denotes  the  rotor  magnetic  axis,  and  the  q  axis  measures  the  angular  displacement 
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of  the  rotor  referenced  from  the  stator  as  winding.  The  q  axis  is  displaced  90°  counter¬ 
clockwise  from  the  d  axis.  The  rotating  magnetic  field  produced  by  the  interaction  of  the 
stator  current  and  the  magnetic  rotor  develops  a  positive  electromagnetic  torque  (TJ, 
which  turns  the  rotor  counter-clockwise  in  the  direction  of  increasing  0^  at  an  electrical 
speed  of  co^-  The  load  torque  (T^)  acts  to  oppose  the  electromagnetic  torque. 

The  following  assumptions  will  be  made  in  our  model  development: 

•  No  demagnetization  of  the  permanent  magnet  due  to  large  currents  or 
temperature. 

•  Eddy  currents  and  harmonic  losses  are  neglected. 

A  loop  analysis  of  the  circuit  shown  in  Figure  3-5  provides  the  motor  voltage 
equations 


Vas=rsL+PK 


(3.2) 


Vfe  =rJt,+pK 


(3.3) 


where  p  is  the  notation  used  for  the  operator  d/dt  and  A,  is  referred  to  as  the  flux  linkage. 
The  sinusoidal  voltages  and  are  90°  phase  shifted.  The  flux  hnkage  is  a  function  of 
the  inductance  of  the  stator  windings  and  the  magnetic  flux  produced  by  the  permanent- 
magnet  rotor.  The  flux  linkage  equations  may  be  expressed  as 

~  ^asaJ'as  "**  ^asbJ’bs  ^asm,  magnetic  rotor  (3.4) 


'k.  =  T.  i 


(3.5) 


and 


{,^bsm  J 


= 


sin  0^ 

[-  cos  0^  J 


(3.6) 


where  is  the  amplitude  of  the  flux  linkages  as  seen  from  the  stator,  0^  is  defined  as 


0,  = 


(3.7) 
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where  co^  is  the  electrical  angular  velocity  and  ^  is  a  dummy  variable  of  integration. 

By  assuming  a  uniform  air  gap  and  identical  windings,  the  mutual  inductance 
between  the  as  and  bs  stator  windings  is  zero  and  the  self-inductances  and  are 
equal  and  denoted  Utilizing  L„  and  equations  (3.4)  through  (3.7),  the  voltage 
equations  (3.2)  and  (3.3)  become 

^as  ^s^as  P(  ^ss^as  (3-8) 

^bs  =  rjbs  +  P(  LJbs  +  (3-9) 

or  in  matrix  form 

V.*.  =  +  pil^Jiaks  +  )  (3. 10) 

The  simulation  of  the  permanent-magnet  machine  is  made  easier  by  transformation 
to  a  reference  frame  fixed  in  the  rotor  by  use  of  the  Pa±’s  transformation  [Ref.  5,  pp.  133- 
160].  In  this  transformation,  the  electrical  frequency  0)^  =  (O^  and  therefore  0g  =  6^  Ibis 
will  enable  the  use  of  a  dc  voltage  input  to  control  the  machine  speed  rather  than  a 
sinusoidal  voltage  from  a  dc-to-ac  inverter.  The  Park's  transformation  for  a  two-phase 
machine  is 

fgsl  _ 

/xj  [sinQr 

or 

(3.12) 
and 

(3.13) 

where  /represents  either  voltage,  current,  or  flux  linkage.  The  subscript  s  denotes  stator 
variables  and  superscript  r  indicates  transformation  to  the  rotor  reference  frame.  By 
substituting  (3.13)  into  (3.10),  the  matrix  voltage  equation  becomes 


sinQr  Jfc 


cosQ^llfbs 


(3.11) 
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(3.14) 


Pre-multiplying  equation  (3.14)  by  and  considering  A,^  as  a  constant,  the  voltage 
equations  may  be  written  in  matrix  form  as 


^  P^ss  ^r^ss  ^qs  ®r^m 

+ 

L  ~®r4.  r,+pL,,\[i^  0 


(3.15) 


Since  0^  =  9^ 


and  equation  (3.15)  becomes 


(3.16) 


^  P^ss  ^r^ss  ir^ 


'(^r^ss  PAyj  JL^'< 


ks  ®Am 
+ 


(3.17) 


The  electromagnetic  torque  is  related  to  the  rotor  speed  by 


T,=(J,  +  J^)-^+  50),  + 


(3.18) 


and  to  motor  current  by 


(3.19) 


Combining  equations  (3.18)  and  (3.19),  we  obtain 


■“  "^m^qs  =  (Jr  +  Jml)  +  ^®r  +  ^motload 


(3.20) 
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where  and  represent  the  inertia  of  the  rotor  and  the  connected  mechanical  load 
respectively,  B  represents  the  damping  coefficient  due  to  the  friction  and  windage  losses, 
and  represents  the  load  torque  delivered  to  the  motor. 

2.  The  Motor  as  a  System  of  State  Variables 

From  equations  (3.17)  and  (3.20),  expressions  for  the  relationship  between  the 
voltage  equations  and  the  electromagnetic  torque  can  be  written  in  terms  of  the  state 
variables  i^^,  ij,  and  co^ 


-0,4-^ct),+— v; 


(3.21) 


dids  ^ 
dt 


+  ^/qs 


(3.22) 


d<s>,  B  ,  fPV  j. 

dt  +  '  \2)  +  '• 


1 

(Jr+Jna) 


T 

^motload 


(3.23) 


By  substituting  the  electrical  time  constant  Xa=LJra,  the  back-emf  constant 


K  -  torque  constant 


into  the  two-phase  system  described 


by  the  linear  differential  equations  (3.21)  through  (3.23),  the  system  becomes  a  three- 
phase  model  represented  by 


(3.24) 


The  constants  k^,  and  are  usually  provided  by  the  manufacturer's  for  a  specific 
motor. 

Two  important  observable  values  are  the  current  drawn  and  the  power  delivered 
by  the  motor.  The  actual  current  drawn  by  the  motor  is  related  to  the  input  by 
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(3.25) 


^dc^dc 


where  and  1^^  are  the  input  voltage  and  current  to  the  motor.  Since  then  the 

current  drawn  by  the  motor  is  given  by 


4c 


(3.26) 


The  actual  power  delivered  by  the  motor  is 


‘  motor 


^motload 


_2> 

PJ 


CO, 


(3.27) 


The  time-domain  block  diagram  for  the  permanent-magnet  dc  machine  represented 
by  equation  (3.24)  and  (3.26)  is  shown  in  Figure  3-6.  The  diagram  represents  the  motor 
component  of  the  vehicle  system  to  be  simulated. 
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Figure  3-6.  Time-domain  block  diagram  implementation  of  the  permanent-magnet 
brushless  dc  motor  given  in  equation  (3.24). 
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IV.  MODELING  THE  VEHICLES  MECHANICAL  CHARACTERISTICS 


A.  THE  FORCES  ACTING  ON  THE  VEHICLE 

The  forces  that  have  been  considered  important  in  modeling  an  electric  vehicle  are 
described  as  the  vehicles  performance  variables  [Ref.  7,  App.  B,  and  Ref.  1,  pp.  46-47]. 
These  performance  variables  are  comprised  of  forces  acting  on  the  vehicle  and  are  used  to 
calculate  the  power  required  to  propel  a  vehicle  at  some  steady  velocity.  The  modeling 
equation  which  describes  this  power  is 

Pmotor  -  power  required  to  drive  the  vehicle 


where,  P/j  =  Power  required  at  the  rear  drive  wheel,  P^  =  Power  required  for  accessories, 
Ej^  =  Efficiency  of  the  accessories,  =  mechanical  and  electrical  efficiency  of  the  drive 
train  and  motor.  By  neglecting  the  accessories,  we  may  express  Pmotor  ^ 


_  Pr  _  v( Fw  +  pR  Phg  Pa) 


(4.2) 


where  v  =  the  vehicle  velocity,  P'w  =  wind  force,  Fg  =  rolling  resistance  force,  F^q  =  hill 
gravitational  force,  and  =  acceleration  force.  These  forces  are  describe  the  total  drag 
on  the  vehicle,  in  Newton's  and  will  be  expanded  upon.  Equation  (4.2)  will  be  referred  to 
as  the  performance  equation. 

1.  The  Wind  Force  (F^) 

This  performance  variable  describes  the  force  of  faction  on  the  vehicle  body  as  it 
moves  through  the  air.  It  is  specified  at  zero  wind  velocity  and  is  a  function  of  the  frontal 
surface  area  of  the  vehicle,  as  well  as  any  protrusions  on  the  vehicle,  air  passages  in  the 
body  and  any  other  components  of  the  vehicle  affected  by  movement  over,  around  or 
through  the  vehicles  shell  or  body.  The  ideal  vehicle  shape  is  a  long  cylinder  with  tapered 
ends.  The  equation  that  describes  the  wind  force  is; 
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where  is  the  air  density  (1.2  kg/m^  was  used  in  all  calculations),  v  is  the  vehicle  velocity 
(m/s),  Q  is  the  coefficient  of  aerodynamic  drag,  and  A  is  the  vehicles  frontal  surface  area 
(m2). 

2.  The  Rolling  Resistance  Force  iJF^ 

The  rolling  resistance  is  primarily  a  function  of  the  tire  friction  on  the  roadway. 
The  coefficient  of  friction  is  directly  related  to  the  type  of  tire  used.  Narrow,  high 
pressure  tires  will  have  a  lower  coefficient  of  friction  than  tires  that  are  wider  and  lower 
pressure.  Bicycle  tires  are  normally  used  in  solar  powered  vehicle  applications.  The 
equation  which  describes  the  rolling  resistance  force  is: 

=  MyC,g  (4.4) 

where  My  =  vehicle  mass  (kg),  =  coefficient  of  rolling  resistance,  and  g  =  gravitational 

constant  (9.8  m/s2). 

3.  The  Gravitational  Force  (F^q) 

The  hill  gravitational  force  affects  the  total  drag  only  when  the  vehicle  is  traversing 
a  hill.  Uphill  equates  to  a  positive  grade  and  downhill  a  negative  grade.  The  equation 
describing  this  force  is: 

FfjQ  =  Myg  sin  0^  (4.5) 

where  0  =  angle  of  the  road  measured  from  the  horizontal  component  of  the  driving  route. 

4.  The  Acceleration  Force  (F^) 

The  acceleration  force  is  also  known  as  the  inertial  force  [Ref.  4,  pp.  20-23].  This 
component  of  the  total  drag  will  have  a  positive  sign  for  acceleration  and  negative  sign  for 
deceleration.  The  acceleration  force  equation  is: 

F^  =  Mya  (4.6) 

where  a  =  vehicle  acceleration  (m/s2).  is  a  dominant  force  in  the  process  of 
regenerative  braking.  This  is  the  process  of  recovering  energy  in  the  form  of  current 
supplied  to  the  battery  as  discussed  in  Chapter  HI. 

5.  The  Load  Torque  to  the  Motor 

The  forces  expressed  in  the  performance  equation  (4.2)  relate  directly  to  the  load 
experienced  by  the  motor.  Realizing  that  a  chain  or  belt  drive  translates  these  forces  to 
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the  motor  tangent  to  the  rotation  of  the  drive  wheel  in  contact  with  the  road,  the  load 
torque  then  becomes  [Ref.  8,  p.  183] 


where  r  is  the  radius  of  the  drive  wheel  and  N  is  the  speed  or  gear  ratio  from  the  drive  axle 
to  the  motor.  Expanding  the  sum  of  vehicle  forces  represented  in  equation  (4.7)  gives 


Figure  4-1  shows  how  the  vehicle  forces  ate  translated  to  the  motor. 


Figure  4-1.  Interaction  of  vehicle  force  with  motor  for  computii^  load  torque  7)^. 

B.  THE  VEHICLE  DRIVE  TRAIN 

The  drive  train  acts  as  the  coupling  mechanism  between  the  motor  and  the  wheels. 
The  most  common  method  used  in  solar  powered  electric  vehicles  are  chain  or  belt  drives. 
The  drive  train  model  will  be  presented  by  assuming  a  belt  drive  configuration.  The  motors 
radial  velocity  in  revolutions  per  minute  is  given  as 
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where  0)^  is  the  rotor  electrical  frequency  in  radians  per  second,  and  P  is  the  number  of 
poles  for  the  specific  motor. 

The  RPM^f  is  converted  to  the  appropriate  linear  velocity  by  multiplication  with 

the  gear  ratio  and  tire  radius.  The  selection  of  the  gear  ratio  depends  on  the  things  such  as 

type  of  terrain  over  which  the  vehicle  is  traveling,  the  desired  vehicle  speed,  the  rated 

speed  of  the  motor,  etc.  For  example,  over  a  level  terrain  a  ratio  of  1:5  may  be  utilized.  If 

on  the  other  hand,  the  terrain  is  uphill,  an  appropriate  gear  ratio  would  be  1:9.  A  higher 

gear  ratio  provides  more  power  for  climbing  hills,  but  also  greatly  reduces  the  maximum 

speed  of  the  vehicle.  Utilizing  equation  (4.9),  the  gear  ratio  G»= — ,  and  the  wheel  radius 

N 

r,  the  vehicle  velocity  is  expressed  by 

Vcarikph)  =  RPM^,  •  •  2 nr  (4.10) 

c^rV  y  mot  ^ 


then  converting  to  miles  per  hour 


Vcar  (  )  =  0.  625 V^ar  (  ) 


(4.11) 


C.  THE  MECHANICAL  COMPONENTS  DESCRIBED  AS  A  SYSTEM 

The  mechanical  components  comprising  the  vehicle  forces  which  translate  into  the 
load  torque  experienced  by  the  motor,  are  presented  in  the  block  diagram  shown  in 
Figure  4-2.  The  input  is  the  angle  of  the  incline  or  decline  in  radians.  If  the  vehicle  is 
traveling  on  level  ground,  then  the  input  would  equal  zero.  The  vehicles  forward 
velocity  component  of  the  performance  equation  (4.2)  is  derived  from  the  drive  train.  The 
velocity  described  in  equation  (4.11)  is  converted  firom  kilometers  per  hour  to  meters  per 
second. 
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The  drive  train  equations  (4,10)  and  (4,11)  are  presented  in  block  diagram  form  in 
Figure  4-3, 


Figure  4-3.  Block  diagram  of  drive  train  conversion  of  motor  rotational  speed  to 

vehicle  forward  velocity. 
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V.  SIMULATION  SOFTWARE  AND  MODEL  DEVELOPMENT 


A.  THE  SIMULATION  SOFTWARE 

Matlab  [Ref.  2]  was  chosen  as  the  simulation  software  for  the  vehicle  system.  The 
Matlab  Toolbox  -  Simulink,  is  the  primary  building  tool  for  designing  and  connecting  the 
dynamic  model  components  of  the  vehicle  simulator.  Simulink  models  both  linear  and 
non-linear  system  components  in  block  diagram  "windows".  These  windows  resemble  the 
block  diagrams  presented  in  many  control  theory  textbooks,  which  allows  the  user  to 
easily  construct  and  modify  a  system  with  a  click  of  a  mouse  button.  The  model 
components  presented  in  the  previous  chapters  are  assembled  in  Simulink  block  diagram 
submodels,  then  connected  to  present  the  overall  vehicle  system.  Figure  5-1  shows  the 
Simulink  model  for  the  entire  solar-powered  electric  vehicle  simulator  system. 

The  user  interface  with  the  vehicle  simulator  is  generated  in  MATLAB  code 
(Appendix  A).  A  user  control  window  initiates  the  program  and  allows  the  user  to  input 
the  vehicles  characteristics  by  clicking  the  mouse  on  a  "pushbutton"  within  the  SPEV 
figure  window.  The  required  input  information  is  displayed  in  the  MATLAB  control 
window.  Once  the  vehicle  data  has  been  entered  the  simulation  may  be  run  and  output 
data  analyzed.  The  specific  vehicle  information  required  for  input  is  discussed  in  Chapter 
VI.  User  instructions  are  provided  in  Appendix  A. 

B.  THE  ELECTRICAL  COMPONENTS  IN  SIMULINK 

In  Chapter  HI,  the  individual  components  which  comprise  the  electrical 
characteristics  of  the  solar  powered  vehicle  were  presented  as  equations  and  block 
diagrams.  These  block  diagrams  will  now  be  presented  in  Simulink  form. 

1.  The  Run  Time  and  Solar  Array  Characteristics 

The  first  components  of  the  vehicle  simulator  are  the  run  time  and  solar  array 
output.  Current  vectors  are  developed  from  the  three  types  of  solar  array  configurations 
presented  in  Chapter  HI,  which  are  the  flat,  ciurved,  and  tracking  array  configurations.  An 
array  type  and  expected  peak  current  output  are  specified  by  the  user,  then  the  associated 
vector  is  read  into  the  Array  look-up  table  block.  The  time  of  day  for  simulation  start  is 
specified  by  the  user  and  can  be  between  0900  and  1700.  Figure  5-2  shows  the  simulink 
implementation  of  the  time  and  array  look-up  table.  The  variable  cloud  cover  block  is 
used  to  simulate  the  loss  of  array  output  due  to  overcast  conditions.  The  user  can  "double 
click"  on  the  block,  opening  it  to 
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Time  in  HoursN —  12:34 


Figure  5-1.  Solar-Powered  Electric  Vehicle  simulator  in  Simulink. 
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of  important  characteristics 


a  slider  gain.  Clicking  on  the  down  or  up  arrow  changes  the  gain  coefficient  which 
produces  changes  in  the  array  current  output  fed  to  the  Electrical  Components  block, 
thereby  simulating  current  degradation  due  to  cloud  cover. 


Cutrentfiom 
Solar  Array 

variable  cloud 
cover 


hours 


Time  in  HoursH — 


12:34 


t 


Run  Time  (seconds) 


Hiis  blodt  cooveits  seconds 
to  hours 


\  This  block  is  a  "Slider  Gidn” 
^ch  mult^lies  die  aoay  current 
a  coefficient  between  0  and  1 


>  Aiiaycunent  to  ffiectricalCom|KHiaib  block 


Figure  5-2.  Simulink  implementation  of  the  Solar  Array. 


2.  The  Electrical  Components  Block 

The  Electrical  Components  block  shown  in  Figure  5-3  takes  input  from  the 
vehicle  throttle,  which  is  a  slider  gain,  current  fed  from  the  array  as  previously  discussed, 
and  current  being  drawn  by  the  drive  motor.  The  output  is  controlled  dc  voltage  to  the 
motor  from  the  battery,  and  a  multiplexed  array  of  outputs  to  be  monitored  in  the  On-Line 
Display.  These  outputs  include: 

•  state-of-charge 

•  total  amps-hours  used 

•  battery  voltage 

•  motor  current  drawn 

•  array  current  delivered 

The  Battery  block  uses  a  look-up  table  to  determine  voltage  output  based  on  state- 
of-charge  as  presented  in  equation  (3.1).  The  user  specifies  the  type  of  battery  to  be  used, 
either  Lead-Acid  or  Silver-Zinc.  The  discharge  curve  for  the  chosen  battery  will  be 
displayed  in  the  Battery  block.  Also  seen  in  Figure  5-3  are  the  contents  of  the  Electrical 
Components  block,  which  will  be  expanded  on. 
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Figure  5-3.  Simulink  representation  of  the  main  electrical  components  of  the  vehicle 

simulator. 
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a.  The  Current  Bus  Block 

The  Current  Bus  block  takes  input  from  the  solar  array  and  the  current 
being  drawn  by  the  motor,  then  uses  these  to  compute  the  battery  state-of  charge  (SOC) 
and  the  amp-hours  (Ahrs)  used.  Figure  5-4  shows  the  expanded  view  of  the  Current  Bus 
block.  As  developed  in  chapter  3,  the  sum  of  the  array  current  and  motor  current  are 
integrated  then,  divided  by  3600  seconds  to  provide  the  amp-hour  usage  output.  Since  the 
total  amp-hour  usage  cannot  be  less  than  zero,  a  switch  is  used  to  "turn  off'  the  amp-hour 
output  at  a  value  of  zero.  The  State-of-Charge  function  block  implements  equation  (3.1) 
by  using  the  amp-hour  input  to  compute  the  battery  state-of-charge  based  on  maximum 
battery  capacity  and  initial  state-of-charge. 

b.  Speed  and  Voltage  Controller  Block 

The  inputs  to  the  Speed  and  Voltage  Controller  block  are  throttle  input 
and  battery  voltage.  The  output  is  controlled  dc  voltage  to  the  motor  to  regulate  speed. 
Figure  5-5  shows  the  expanded  Simulink  block.  The  throttle  input,  from  a  slider  gain  is 
sent  to  the  voltage  ratio  function  block  which  outputs  a  coefficient  from  0  to  1.  The 
throttle  coefficient  is  rate  limited  before  being  multiplied  by  the  battery  voltage.  The 
mayimiim  voltage  to  the  motor  is  set  by  a  saturation  block  at  1(X)  volts.  The  rate  liimted 
or  "throttled"  voltage  is  then  sent  to  the  motor. 

3.  The  Drive  Motor  Block 

The  inputs  to  the  Drive  Motor  block  are  controlled  voltage  from  the  Electrical 
Components  block  and  vehicle  load  The  vehicle  load  will  be  discussed  in  following 

sections.  The  outputs  are  RPM  to  the  drive  train  and  a  multiplexed  array  of  RPM,  motor 
load,  efficiency  and  current  drawn  by  the  motor  to  the  On-Line  Display.  The  Drive  Motor 
block  is  presented  in  Figure  5-6  and  its  contents  are  shown  in  Figure  5-7. 

a.  The  System  Efficiency  Block 

The  System  Efficiency  block  consists  of  a  two  dimensional  look-up  table 
with  inputs  of  motor  rpm  and  vehicle  load,  and  switches  for  minimum  and  maximum 
efficiency.  The  efficiency  limits  are  set  by  the  minimiun  and  maximum  efficiency  for  the 
specific  motor  being  used.  The  System  Efficiency  block  is  presented  in  Figure  5-8.  The 
actual  vehicle  load  is  divided  by  the  efficiency  output  in  the  Load/Eff  function  block.  The 
output  of  this  fimction  block  relates  to  the  T^^^  introduced  in  equation  (3.18),  and  used 
in  equation  (3.27)  to  determine  motor  power  out.  This  is  the  vehicle  load  torque  seen  by 
the  motor. 
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Figure  5-4.  Simulink  representation  of  Ahrs  and  SOC  output  from  current  input 
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Figure  5-5.  Simulink  representation  of  Speed  and  Voltage  control. 
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Figure  5-6.  Simulink  diagram  of  the  Drive  Motor  block. 
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b.  The  PMBDC  Block 

The  PMBDC  block  represents  the  permanent-magnet  brushless  dc  motor 
used  in  the  vehicle  and  implements  equation  (3.24)  in  Simulink  form.  The  Simulink 
diagram  of  the  motor  is  presented  in  Figure  5-9.  The  magnitude  of  the  dc  voltage  through 
input  2  directly  controls  the  desired  rpm  of  the  motor.  The  actual  motor  speed  is  of 
course  affected  by  the  load  torque  introduced  by  the  mechanical  components  of  the 
vehicle  and  the  inertial  load  The  speed  in  rpm  is  output  to  the  Drive  Train  block  of 
the  simulator.  The  motor  current  drawn,  rpm  and  load  torque  are  sent  to  the  On-Line 
display  to  be  monitored  during  simulation. 

C.  THE  MECHANICAL  COMPONENTS  IN  SIMULINK 

In  Chapter  IV,  the  individual  components  which  incorporate  the  mechanical 
characteristics  of  the  solar  powered  vehicle  were  presented  as  equations  and  submodels. 
The  Simulink  form  of  these  submodels  will  now  be  presented  and  their  operation 
explained. 

1.  The  Vehicle  Drive  Train  Block 

The  Vehicle  Drive  Train  block  uses  the  rpm  input  to  convert  rotational  velocity  of 
the  motor  to  the  linear  velocity  of  the  vehicle.  The  Simulink  system  through  which  this  is 
accomplished  is  presented  in  Figure  5-10.  The  rpm  delivered  from  the  Drive  Motor  block 
through  input  1  is  divided  by  a  specified  gear  ratio  and  mechanical  efficiency,  then 
multiplied  by  the  circumference  of  the  drive  wheel.  This  provides  a  linear  velocity  in 
meters  per  minute.  The  linear  velocity  is  converted  to  meters  per  second,  for  use  in  the 
Vehicle  Characteristics  block  in  determination  of  forces,  and  miles  per  hour.  The  miles 
per  hour  variable  is  integrated  to  provide  the  total  miles  traveled.  The  speed,  in  miles  per 
hour,  and  total  distance  traveled  are  sent  to  the  On-Line  Display  and  monitored  as  the 
speedometer  and  odometer.  The  total  miles  traveled  is  also  sent  to  the  Vehicle 
C/iaractemftcs  block  for  force  computations. 

2.  The  Vehicle  Characteristics  (Load  Torque)  Block 

The  Vehicle  Characteristics  (Load  Torque)  block  is  simply  the  implementation  of 
the  performance  equation  (4.2)  and  the  associated  torque  equations  detailed  in  section  4.5. 
The  inputs  to  the  Vehicle  Characteristics  block  are  the  velocity  in  meters  per  second,  total 
miles  driven,  and  the  voltage  delivered  to  the  motor.  Each  of  the  forces  developed  by  the 
motion  of  the  vehicle  are  summed  and  sent  to  the  switch  block.  At  this  switch  block  the 
total  force  is  either  zero,  if  no  voltage  and  no  velocity  are  present,  or  the  sum  of  all  forces 
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Figure  5-9.  Simulink  representation  of  a  Permanent-Magnet  Brushless  DC  Motor 

(PMBDC). 
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present.  Then,  since  the  total  force  is  actually  applied  at  the  point  where  the  wheel 
contacts  the  road,  the  force  is  multiplied  by  the  wheel  radius  and  divided  by  the  gear  ratio, 
providing  the  vehicle  load  torque  seen  by  the  motor.  The  complete  contents  of  the 
Vehicle  Characteristics  block  are  shown  in  Figure  5-11.  Each  of  the  force  components 
within  the  vehicle  characteristics  block  will  be  presented  and  expanded  upon. 

a.  The  Rolling  Resistance  Block 

The  contents  of  the  Rolling  Resistance  block  are  shoAvn  in  Figure  5-12. 
The  input  to  the  block  is  the  linear  velocity  of  the  vehicle  and  the  hill  force.  The  switches 
shown  have  three  purposes:  1)  if  the  vehicle  velocity  is  zero,  then  the  rolling  resistance 
force  is  zero,  2)  if  the  vehicle  velocity  is  greater  than  zero,  the  rolling  resistance  force  is 
applied  against  the  direction  of  motion,  and  3)  if  the  hill  force  is  greater  than  the  rolling 
resistance  force,  and  the  vehicle  is  moving  backwards,  the  rolling  resistance  force  is 
applied  in  the  opposite  direction  of  motion.  The  rolling  resistance  force  is  a  constant  and 
not  a  function  of  the  vehicle  velocity. 

b.  The  Hill  Gravitational  Block 

The  Hill  Gravitational  block  computes  the  force  on  the  vehicle  as  a 
function  of  hill  angle.  If  the  vehicle  is  traveling  up-hill,  the  force  is  positive  and  increases 
the  total  load  on  the  motor.  If  the  vehicle  travel  is  down-hill,  the  force  is  negative  and 
consequently  decreases  the  total  load  on  the  motor.  If  the  vehicle  is  driving  over  level 
terrain,  then  the  hill  force  is  zero.  This  component  plays  a  significant  role  in  the 
regenerative  action  of  the  motor.  If  the  hill  angle  is  large  enough  to  overcome  the  other 
forces  applied  to  the  motor,  then  the  load  of  the  vehicle  will  "drive"  the  motor,  which  in 
turn  win  "generate"  current  for  battery  charging.  Computation  of  the  hiU  force  in 
Simulink  is  shown  in  Figure  5-13. 

c.  The  Acceleration  Force  Block 

The  input  to  the  Acceleration  Force  block  is  the  linear  velocity  of  the 
vehicle.  The  velocity  is  differentiated  to  obtain  the  vehicle  acceleration  used  in 
determining  the  accelerating  or  decelerating  force.  As  with  the  hUl  force,  the  acceleration 
force  also  plays  a  significant  role  in  the  regenerative  action  of  the  motor.  If  the  voltage 
applied  to  the  motor  is  decreased,  the  vehicle  will  begin  to  slow  down.  This  action  causes 
a  large  deceleration,  or  negative  force,  to  be  applied  to  the  motor,  which  in  turn  creates  a 
significant  regeneration  of  current  to  the  battery  for  charging. 
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Figure  5-12.  Simulink  representation  of  the  Rolling  Resistance  force. 
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This  force  is  only  applied  during  transient  operation  and  has  no  steady-state  component. 
Figure  5-14  shows  the  Simulink  representation  of  the  acceleration  force  computation. 


Vehicle 

Mass 


Figure  5-14.  Simulink  representation  of  the  Acceleration  Force  computation. 
d.  The  Wind  Force  Block 

Figure  5-15  shows  the  Simulink  representation  of  the  wind  force.  The 
input  to  the  Wind  Force  block  is  the  linear  velocity  of  the  vehicle. 
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Figure  5-15.  Simulink  representation  of  the  Wind  Resistance  computation. 


This  force  is  directly  proportional  to  the  square  of  the  vehicles  linear  motion.  The  faster 
the  vehicle  goes,  the  greater  the  force. 


40 


D.  THE  VEHICLE  ON-LINE  DISPLAY  IN  SIMULINK 

The  On-Line  Display  block  takes  inputs  from  various  components  of  the 
simulation  program.  These  inputs  are  all  feed  into  a  simulink  multiplex  block  and  then 
sent  to  a  Madab  s-function  called  spevanl.m.  The  function  spevan2.m  is  included  in 
Appendix  B  with  the  other  Matlab  functions  used  in  the  vehicle  simulation  program.  The 
On-Line  Display  block  and  its  components  is  presented  in  Figure  5-16.  All  inputs  are 
clearly  labeled. 

Once  the  simulation  routine  has  been  initiated  by  the  user,  the  On-Line  Display 
window  will  appear  at  the  bottom  of  the  monitor.  The  outputs  to  the  screen  are  clearly 
identified.  Upon  completion  of  a  simulation  run  a  "double-click"  on  the  Plots  block  will 
provide  a  profile  of  important  data  and  average  values  of  certain  characteristics.  These 
average  value  components  will  be  detailed  in  the  following  chapter. 
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Figure  5-16.  On-Line  Display  block  and  contents  in  Simulink. 
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VI.  TESTING  AND  EVALUATION  OF  THE  SYSTEM  MODEL 


A.  OVERVIEW  OF  THE  TESTING  APPROACH 

The  primary  question  to  ask  in  testing  the  vehicle  simulation  program  is,  "Does  the 
power  out  of  the  motor  equal  the  power  required  to  overcome  the  forces  acting  on  the 
vehicle  for  a  specific  speed?"  This  chapter  will  be  concerned  with  answering  this  question 
over  three  types  of  terrain  profiles  -  flat,  single  hill,  and  variable  elevation. 

I.  Specification  of  the  Initial  Simulation  Parameters 

The  vehicle  electrical  and  mechanical  simulation  parameters  where  compiled  from 
actual  solar  powered  vehicles  built  by  teams  entered  in  the  1990  and  1993  World  Solar 
Challenge. 

The  mechanical  characteristics  were  complied  from  a  vehicle  built  and  raced  by  the 
CALSOL  solar  powered  vehicle  race  team  from  the  University  of  Catifomia  at  Berkeley. 
Table  6-1  shows  the  performance  characteristics  used  for  the  vehicle  simulator  testing. 
These  characteristics  also  serve  as  the  default  values  for  initiation  of  all  simulations. 


Coefficient  of  Drag  (Cn) 

.16 

Frontal  Area  (A) 

1.2  m2 

Vehicle  Mass  (Mv) 

295.8  kg 

Coefficient  of  Rolling  Resistance  (Cr) 

.0105 

Rear  Wheel  Radius  (r) 

.26  m 

Gear  Ratio  (1:N) 

1:5 

Drive  Train  Efficiency 

98% 

Table  6-1.  Performance  characteristics  used  in  the  vehicle  simulation. 

A  permanent-magnet  brushless  dc  motor  built  by  Unique  Mobility,  Inc.,  is  used  in 
this  vehicle  simulation.  The  characteristics  of  the  motor  used  in  the  simulation  are 
presented  in  Table  6-2.  Figure  6-1  shows  the  efficiency  map  provided  for  the  motor.  This 
efficiency  map  is  loaded  into  the  program  upon  initialization  of  the  motor  parameters  for 
use  during  simulation.  As  with  the  vehicle  characteristics,  the  motor  parameters  used  for 
simulation  testing  also  serve  as  the  default  values  for  simulator  initialization. 
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Unique  Mobility  DR127s/CR10-100 

Voltage  In 

50  volts 

Rated  Power 

Max  No-Load  Speed 

Number  of  Poles  (P) 

24 

Damping  Coefficient  (B) 

0 

Rotor  Inertia  (J) 

Electrical  Time-Constant  (T^) 

6.944e-4 

Winding  Resistance  (ra) 

3.6e-2  ohms 

Back-EMF  Constant  (k^) 

20  V/Krpm 

Torque  Constant  (k^) 

.191Nm/A 

Table  6-2.  PMBDC  Motor  parameters  used  for  simulation  testing. 


Efficiency  plot  for  Unique  DR1 27/1 0 


Figure  6-1.  Motor  efficiency  plot 
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The  battery  chosen  for  simulator  testing  is  a  Lead- Acid  type,  with  a  capacity  of  35 
Ahrs  and  initial  state-of-charge  of  100%.  A  tracking  solar  array  will  be  used,  with  a  peak 
current  output  of  11  amps.  This  battery  and  array  configuration  are  the  default 
components  used  in  simulation  initialization. 

B.  MODEL  SIMULATION 

The  three  types  of  terrain  considered  during  testing  are  flat,  single  hill,  and 
variable  terrain  which  contains  numerous  positive  and  negative  grades.  The  following 
outputs  will  be  presented: 

•  driving  terrain  profile 

•  vehicle  speed 

•  sum  of  the  forces  on  the  vehicle 

•  current  drawn  by  the  motor 

•  torque  seen  by  the  motor 

•  power  delivered  by  the  motor 

•  power  required  at  the  rear  drive  wheel 
1.  Simulation  Over  Flat  Terrain 

The  first  driving  profile  for  simulation  is  over  flat  terrain  as  shown  in  Figure  6-2. 
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Figure  6-2.  Plot  of  flat  terrain  driving  profile. 

The  simulator  system  was  initialized  with  the  parameters  specified  in  the  previous  section. 
Cloud  cover  was  neglected,  so  peak  current  from  the  solar  array  was  1 1  amps. 
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The  first  outputs  examined  are  the  vehicle  speed  and  total  force  on  the  vehicle. 
The  theoretical  speed  is  calculated  from  a  quadratic  equation  based  on  the  relationship 
between  the  energy  delivered  by  the  motor  and  the  energy  required  to  overcome  the 
vehicle  forces.  This  is  expressed  as 


A  +  B( rpm^t )  =  C  -  D( rpm^, ) 

_  MygCff  _  p^CjyA2n^r^ 

NE^  ’  60^N^E^ 


(6.1) 


With  data  provided  in  Tables  6-1  and  6-2,  we  get 


+ 


0. 1055 


0.  000000607 


^  (6.2-  263.9)  _  ^ 
0.  000000697 


Then  using  the  quadratic  equation 


(6.2) 


rpm^ot 


_  -b  ±  -  4ac 

2a 

=  -153770  or  2404.  5 


(6.3) 


The  motor  rpm  in  (6.3)  equates  to  a  vehicle  speed  of  =  29.4  mph  or  13.1  m/s.  lliis 
value  is  within  95%  of  the  simulated  value  presented  in  Figure  6-3.  The  slight  variation  is 
most  likely  caused  from  not  considering  the  motor  efficiency  in  the  calculation. 

Using  a  speed  of  13.1  m/s,  the  theoretical  total  steady  state  force  acting  on  the 
vehicle  is  derived  from  equation  (4.2)  as 


Pvehick  =  K  Pr 

=  MygCji  -t-  —  p^CjjAv^ 

=  (295.8)(9.  8)(0.  0105)  +  -(1.  2)(0. 16)(1.  2)(13. 1^ ) 

2 

=  50.  2  (  Newtons ) 

(6.4) 
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Figure  6-3.  Plots  of  the  steady  state  speed  and  total  force  acting  on  the  vehicle  for 

flat  terrain  driving 
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The  acceleration  force  (F^)  and  hill  gravitational  force  (F^^)  are  zero  for  flat  steady  state 
driving.  The  force  given  by  (6.4)  is  within  96%  of  the  simulated  steady-state  value  of  the 
total  force  displayed  in  Figure  6-3. 

Also  observed  during  the  simulation  are  the  current  drawn  by  the  motor  and  the 
torque  delivered  to  the  motor.  Using  equation  (4.8),  the  theoretical  force  provided  in 
(6.4),  and  data  given  in  Table  6-1,  the  theoretical  load  torque  delivered  to  the  motor  is 


{0.  26)(50,2) 
(5)(0.98) 


=  2.6 ( Newton  -  meters) 


(6.5) 


Utilizing  the  motor  rpm  from  (6.3)  and  the  load  torque  from  (6.5),  the  motor  efficiency 
may  be  determined  from  the  efficiency  plot  in  Figiue  6-1.  The  efficiency  is  approximately 
75%.  Dividing  (6.5)  by  the  efficiency,  the  load  torque  delivered  to  the  motor  is 

Tn^dnnA  =  — =  S .  47 ( NewtoH  -  meters)  (6.6) 


Then  utilizing  (6.6)  and  equation  (3.24)  the  theoretical  current  drawn  by  the  motor  is 


As  seen  in  Figure  6-4,  the  theoretical  load  torque  and  motor  current  drawn  correspond 
closely  to  the  steady-state  output  from  the  simulated  model.  Both  values  are  within  97% 
of  expected  results. 

The  most  important  comparison  to  make  is  the  power  required  to  drive  the  vehicle 
versus  the  power  delivered  by  the  motor.  The  theoretical  power  required  is  computed 
from  the  performance  equation  (4.2).  Utilizing  the  vehicle  velocity  from  (6.3)  and  the 
total  force  on  the  vehicle  given  in  (6.4),  the  power  required  at  the  rear  wheels  is 

P^=v(F^  +  Fj,)  =  (13. 1)(50. 2)  =  657. 62(  Watts)  (6.8) 


Dividing  (6.8)  by  the  motor  efficiency  and  the  drive  train  efficiency,  the  power 
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Figure  6-4.  Plot  of  the  motor  current  and  the  load  torque  delivered  to  the  motor  for 

flat  terrain  driving. 
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required  from  the  motor  is 


R _  _ 


657.  62 
(0.75)(0.98) 


=  894.  72  W 


The  output  plots  of  the  power  delivered  by  the  motor  and  the  power  required  at 
the  rear  wheels  with  consideration  for  system  efficiencies  are  presented  in  Figure  6-5.  The 
simulated  steady-state  values  are  within  94%  of  the  computed  power  in  (6.9). 

2.  Simulation  Driving  Over  a  Single  Hill  Profile 

The  next  simulation  was  up  a  single  hiU.  Again,  the  array  current  was  a  constant 
1 1  amps.  The  terrain  profile  for  this  simulation  is  shown  in  Figure  6-6.  Using  (6.1),  with 


^  ^  MygCg  ^  MygsinQ^ii 


(6.10) 


the  speed  for  the  terrain  profile  with  0^;  =  0.  03789  radians  is 


2  0. 1055  (28.62-  263.9)  _ 

+  Q  000000697  0.  000000697 


(6.11) 


(6.12) 


The  motor  rpm  in  (6.12)  equates  to  a  vehicle  speed  of  =  2<5,9  mph  or  12.0  m/s.  This 
value  is  within  94%  of  the  simulated  value  presented  in  Figure  6-7. 

Using  a  speed  of  12.0  m/s,  the  theoretical  total  steady  state  force  acting  on  the 
vehicle  is 

Vehicle  ~  ^R  ^HG  ~  ^vS^'R  +  ^v8  ^hUl 

=  (295.  8)(9.8 )[(0.  0105)  +  sin(0.  03789)]  +  |  f  2)( 0. 16 )( 1.  2)(12.  0^  ) 
=  156.  8 ( Newtons) 

(6.13) 
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Figure  6-5.  Plots  of  the  power  required  at  the  rear  wheel  and  power  delivered  hy 

the  motor  for  flat  terrain  driving. 
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Single  Hill 


Figure  6-6.  Plot  of  single  hill  terrain  profile . 
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Figure  6-7.  Plots  of  the  vehicle  speed  and  the  total  force  acting  on  the  vehicle  for  a 

single  hill  terrain  profile. 
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The  simulated  force  in  Figure  6-7  is  within  97%  of  the  computed  force  in  (6.13). 

Using  equation  (4.8),  the  force  from  (6.13),  and  data  from  Table  6-1,  the 
theoretical  torque  delivered  to  the  motor  for  this  terrain  profile  is 

T  =  +  Fr  +  Fhg) 

NE^ 

^  (0.  26)(156.8) 

(5)(0.98) 

=  8.  32 ( Newtons  -  meters ) 


I  Trill  zing  (6.12),  (6.14),  and  the  efficiency  plot  in  Figure  6-1,  the  motor  efficiency  is 
approximately  85  %.  By  dividing  (6.14)  by  the  motor  efficiency,  the  theoretical  load 
torque  delivered  to  the  motor  is 

Tnu>tioad  =  7^  =  9.78 ( Newton  -  meters)  (6.15) 

U,  oj 


Then  utilizing  (6.15)  and  equation  (3.24)  the  theoretical  current  drawn  by  the  motor  is 


51.  2(  Amps  ) 


(6.16) 


As  seen  in  Figure  6-8,  the  theoretical  load  torque  and  current  are  within  92%  of  the  values 
obtained  from  the  simulation. 

For  the  power  comparison,  the  theoretical  power  required  for  this  driving  profile  is 
computed  from  the  performance  equation  (4.2).  Utilizing  the  total  force  on  the  vehicle 
given  in  (6.13)  and  the  vehicle  velocity,  the  power  required  at  the  rear  wheels  is 

Pg=v(F„  +  Fg  +  Fg^)  =  (12.0)( 156.8)  =  1881.6(Watts)  (6.17) 

Dividing  (6.17)  by  the  motor  efficiency  and  drive  train  efficiency,  the  power  required  from 
the  motor  is 


_ 


1881.6 

(0.86)(0.98) 


=  2258. 5(  Watts) 


(6.18) 


54 


current  drawn  by  motor  (A) 


Load  torque  in  (Nm) 


miles  driven 


Figure  6-8.  Plots  of  the  current  drawn  by  the  motor  and  the  load  torque  delivered 
to  the  motor  for  a  single  hill  driving  profile. 
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The  output  plots  of  the  power  delivered  by  the  motor  and  the  power  required  at  the  rear 
wheels  with  consideration  for  system  efficiencies  are  presented  in  Figure  6-9.  The 
simulated  values  are  within  92%  of  the  expected  value  of  power  computed  in  (6.18).  If 
the  vehicle  speed  had  been  computed  with  consideration  of  motor  efficiency,  there  would 
have  been  less  variation  between  theoretical  and  simulated  results. 

3.  Simulation  Driving  Over  Variable  Elevation  Terrain 

The  next  simulation  is  a  variable  or  "hiUy"  terrain.  Solar  array  current  was  a 
constant  1 1  amps.  The  terrain  profile  for  this  simulation  is  shown  in  Figure  6-10. 

Since  the  simulation  has  been  shown  to  correspond  closely  with  expected  values 
for  the  flat  and  single  hill  driving  profiles,  the  plots  for  this  variable  terrain  will  be 
presented  without  theoretical  comparison.  As  seen  in  Figures  6-11  through  6-13,  the 
outputs  are  consistent  with  the  results  comparisons  made  in  the  two  previous  driving 
profiles.  Observation  of  Figure  6-13  shows  that  the  power  delivered  equals  the  power 
required.  This  would  indicate  that  the  simulation  model  does  accurately  portray  an 
electric  vehicle  driving  over  a  variable  terrain. 

Also  observed  in  Figure  6-13  is  the  regenerative  action  of  the  motor  when  the 
forces  acting  against  the  vehicle  motion  are  overcome  by  forces  which  act  to  propel  the 
vehicle  forward.  This  creates  a  current  regeneration  action  from  the  motor  which  is  used 
by  the  system  to  charge  the  battery,  increasing  the  battery  state-of-charge. 

4.  The  Optimum  Balance  Between  Motor  Current  and  Speed 

Since  the  objective  of  this  simulation  model  is  to  optimize  the  design  of  the 
vehicle,  the  next  test  should  be  to  find  the  most  efficient  relationship  between  current 
drawn  by  the  motor  and  speed  of  the  vehicle.  This  test  was  driven  over  flat  terrain  and  the 
voltage  to  the  motor  was  stepped  in  increments  of  10  volts  using  the  throttle  control.  At 
each  voltage,  the  system  was  allowed  to  settle  to  steady  state,  then  the  current  being 
drawn  by  the  motor  and  the  vehicle  speed  were  recorded.  Table  6-3  shows  the  data 
accumulated  from  this  experiment 

Since  the  current  drawn  by  the  motor  is  a  function  of  the  motor  efficiency,  there 
should  be  a  range  of  speed  in  which  the  motor  is  most  efficiently  using  the  energy 
supplied.  By  plotting  the  current  as  a  function  of  throttle  voltage  and  speed,  the  most 
efficient  operating  points  may  be  easily  identified.  By  examining  the  output  relationship 
shown  in  Figure  6-14,  the  optimum  driving  speed  for  this  particular  vehicle  configuration 
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Power  required  at  the  rear  wheel 
considering  motor  efficiency  (W) 


miles  driven 


Power  delivered  to  the  rear  wheels  (W) 
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Figure  6-9.  Plots  of  the  power  required  at  the  rear  wheel  and  the  power  delivered 
by  the  motor  for  single  hill  driving  profile 
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Figure  6-11.  Plots  of  the  vehicle  speed  and  the  total  force  acting  on  the  vehicle  for  a 

variable  terrain  profile. 


current  drawn  by  motor  (A) 


Load  torque  (Nm) 


Figure  6-12,  Plots  of  the  current  drawn  by  the  motor  and  the  load  torque  delivered 
to  the  motor  for  a  variable  terrain  profile. 
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power  power 


Power  required  at  the  rear  wheel 
considering  motor  efficiency  (W) 


Power  delivered  to  the  rear  wheel  (W) 


Figure  6-13.  Plots  of  (he  power  required  at  (he  rear  wheel  and  the  power  delivered 
by  the  motor  for  a  variable  terrain  profile. 
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would  be  between  25  and  30  mph.  This  test  could  also  be  used  to  determine  the  optimum 
driving  speeds  for  variable  terrain's  and  different  gear  ratios. 


Throttled  Voltage  (V) 

Motor  Current(A) 

Vehicle  Speed  (mph) 

14 

26.5 

7.5 

24 

21.4 

13.5 

34.2 

18.9 

19 

44.2 

17.4 

24.9 

54.3 

17.35 

30.2 

64.4 

18.5 

35.5 

74.5 

20.3 

40.1 

84.5 

22.6 

44.9 

94.7 

24.9 

49.1 

100 

26.1 

51.2 

Table  6-3.  Motor  Current  and  Vehicle  Speed  for  Incremental  Voltages. 

Efficient  vehicle  configurations  can  be  achieved  by  changing  the  individual  vehicle 
parameters  necessary  to  optimize  the  current  -  speed  relationship. 
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vn.  CONCLUSIONS 


There  is  a  need  for  efficient  analysis  of  alternative  energy  propulsion  systems  for 
the  automobiles  of  tomorrow.  By  developing  systems  that  quickly  analyze  a  particular 
vehicle  configuration,  the  system  designer  can  save  the  time  and  money  required  to 
actually  build  full  scale  models.  This  thesis  has  provided  a  method  for  the  efQcient 
evaluation  of  electricity  as  an  alternative  energy  source  in  automobile  propulsion. 

A.  CURRENT  APPLICATION 

The  original  goal  of  this  thesis  was  to  model  a  solar  powered  electric  vehicle  for 
use  in  design  analysis  of  vehicles  configured  under  the  rules  of  the  World  Solar  Challenge. 
That  task  has  been  successfully  completed. 

Presentation  of  this  thesis  was  made  to  the  Solar  Phantom  Team  at  Rose-Hulman 
Institute  of  Technology  in  Terra  Huete,  Indiana  on  March  6,  1995.  The  system  evaluation 
and  output  presented  in  the  previous  chapters  were  verified  by  observations  made  by  the 
Rose-Hulman  team  during  previous  solar  powered  races  in  the  Solar  Phantom  I  and  11. 
The  team  is  presently  constructing  the  Solar  Phantom  III,  which  has  a  radically  different 
design  than  the  previous  two  cars.  The  simulation  model  developed  in  this  thesis  wiU  be 
used  to  develop  a  daily  driving  strategy  for  the  Solar  Phantom  Team  when  they  participate 
in  Sunrayce  95  on  June  20,  1995.  This  strategy  should  assist  the  team  in  efficiently  using 
the  energy  available  from  the  solar  array,  battery  and  motor  at  variable  speeds  and  over 
variable  terrain. 

B.  FUTURE  WORK 

This  system  is  by  no  means  complete.  Future  work  could  include  improvement  in 
the  areas  of: 

•  modeling  temperature  effects  on  the  motor  and  solar  cells 

•  modeling  the  motor  to  account  for  harmonic  losses 

•  modeling  a  friction  braking  mechanism 

•  adding  regenerative  current  control 

•  incorporating  the  effects  of  auxiliary  components  in  power  consumption 

•  incorporate  a  random  number  generator  for  the  cloud  cover  coefficient  to 
simulate  variability  of  sunlight  to  the  solar  array 

Other  work  which  could  make  this  system  more  ^plicable  to  industry  would  be  to 
develop  the  models  for  major  components  such  as  the  battery,  motor,  vehicle 
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characteristics,  and  solar  array  as  individual  "plants"  and  evaluate  them  with  a  digital 
Signal  Processing  and  Control  Engineering  (dSPACE)  system.  The  plant  models  could  be 
built  in  actual  hardware,  then  connected  to  the  simulation  program  through  the  d-Space 
interface.  Through  this  connection,  the  simulation  program  may  be  used  to  operate  the 
hardware  plant.  The  systems  designer  may  observe  how  hardware  being  considered  for 
use  in  the  vehicle  design,  reacts  under  actual  operating  conditions.  The  d-Space  model 
would  be  much  more  cost  efficient  for  actual  component  evaluation  than  construction  of 
full  scale  models. 
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APPENDIX  A.  SIMULATOR  USER  INSTRUCTIONS 


1.  The  Solar/Electric  Vehicle  Simulator  will  only  run  from  a  Madab  control  window.  The 
specific  version  of  Madab  is  at  least  4.1  with  Simulink  version  1.3c.  Ensure  that  these  are 
loaded  on  your  system  before  attempting  to  run  the  simulator. 

2.  Once  in  the  Madab  control  window,  type  -  solar.  The  user  control  window  and  the 
default  terrain  profile  figures  should  appear  in  the  upper  right  comer  of  the  monitor. 
Figure  A-1  shows  the  user  control  window  for  the  simulator. 


-  I _ Solar/Eleclric  Vehicle  Simulator 

FQe  Edit  \riiidow  Help 


1 

Figure  A-1.  SPEV  simulation  window 


Figure  A-2  shows  the  defaidt  driving  terrain  profile  loaded  upon  simulator  system 
initialization.  Each  button  in  the  user  control  window  invokes  a  specific  function.  Theses 
functions  are  explained  below 

a)  Load  finitial  Data:  Loads  the  initial  default  values  for  all  electrical  and 
mechanical  parameters.  Once  loaded,  aU  values  are  displayed  in  the  Madab  control 
window. 

b)  Run  Time:  Modifies  the  total  simulation  time,  time  of  day  to  start  driving,  and 
integration  time  step. 

c)  Motor  Parameters:  Allows  the  user  to  modify  or  change  the  motor 
characteristics,  or  implement  a  new  motor.  Input  the  requested  data  at  the  prompts  in  the 
Madab  control  window. 
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d)  Vehicle  Parameters:  Allows  any  of  the  vehicles  mechanical  characteristics  to 
be  modified.  Input  requested  data  at  the  prompts  in  the  Madab  control  window. 


-  Tenain  Profile 

File  Edit  ^^iidow  He)^ 


Figure  A-2.  Default  driving  terrain  profile 


e)  Array  Type:  Change  between  tracking,  flat,  or  curved  solar  array 
configurations  and  specify  the  peak  array  current 

f)  Battery  Type:  Change  between  Lead-Acid  or  Silver-Zinc  battery  packs. 

g)  Terrain  Profile:  Allows  the  user  to  choose  one  of  three  terrain  profiles.  At 
the  prompt  in  the  Madab  control  window  choose  flat  single  hiU,  or  load  a  file  of  variable 
terrain  data.  If  the  variable  terrain  is  chosen,  the  data  must  be  entered  from  a 
"filename.dat"  text  file  with  two  columns  of  information.  Column  one  has  distance  in 
increments  of  .1  miles  and  column  two  contains  the  corresponding  elevation  in  feet  above 
sea  level  Each  time  a  terrain  profile  is  changed,  a  modified  terrain  profile  window  is 
displayed. 

h)  Show  All  Data:  Provides  an  updated  display  of  all  important  parameters 
loaded  into  the  simulation  program.  This  function  useful  in  verifying  changes  in  system 
parameters. 

i)  Motor  Efficiency  Plot:  Provides  a  3-dimensional  plot  of  the  motor  efficiency 
as  provided  by  the  manufacturer.  The  user  must  supply  the  efficiency  matrix  for  the 
specified  motor.  The  columns  correspond  to  increasing  load  torque,  while  the  rows 
correspond  to  increasing  motor  rpm. 

j)  Welcome:  Shows  a  brief  welcome  to  the  simulator. 
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k)  Help:  Provides  a  brief  help  on  running  a  simulation  from  the  simulation 
window. 

l)  Close:  Closes  the  user  window. 

m)  Open  Simulation  Window;  Opens  the  primary  vehicle  simulation  window. 
Once  all  previous  parameter  input  requirements  have  been  meet,  simply  click  on  the 
"Simulation  -  Start"  button. 

3.  Once  the  simulation  has  begun,  the  data  output  window  shown  in  figure  A-3  will  be 
displayed  at  the  bottom  left  of  the  monitor. 


— 

Output 

Fi] 

le  Edit 

Wmdow  Help 

Motor  Current  (amps): 

• 

Aim  Used: 

Array  Current  (amps):  TtPM: 

Run  time  (secs):  Hill  Grade  (%): 

Odometer  (miles): 

Speedometer  (mph): 

System  Efficiency  (%): 

Mph 

^  30 

SOC(%); 

Battery  Vohage  (vohs): 
Vehicle  Load  (N-m): 

\ 

45 

Figure  A-3.  Simulation  On-Line  Display  window 


The  output  indicators  are  self-explanatory  and  wiU  provide  continuously  updated  data 
display  throughout  the  simulation. 

4.  The  simulation  may  be  stopped  at  any  time  simply  by  clicking  on  the  "Simulation  - 
Stop"  button  in  the  main  Simulink  simulation  window. 

5.  By  double-clicking  on  the  "Plots"  box  in  the  Simulink  simulation  window,  plots  of  the 
important  vehicle  characteristics,  as  well  as  average  values  of  output  data  wiU  be  displayed 
in  four  separate  figure  windows.  Each  of  these  windows  may  be  closed  by  clicking  on  the 
"File  -  Close"  indicator. 
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APPENDIX  B.  PROGRAM  FUNCTIONS 


A.  FUNCTION  SOLAR.M  FOR  MAIN  PROGRAM 

%  Solar  Powered  Electric  Vehicle  Simulator 
%  Masters  Thesis  EE  Naval-Post  Graduate  School,  Monterey 
%  This  system  will  simulate  a  solar/electric  powered  vehicle 
%  for  any  type  of  vehicle  mechanical  or  electric  characteristics. 
% 

%  This  system  may  be  used  to  optimize  vehicle  configurations. 
%  each  function  is  explained  separately... 

% 


%  This  function  initiates  the  system  user  control  window  and 
%  the  default  terrain  profile. 

% 


%  This  system  was  developed  by  Steven  Roerig,  LT,  USN 
%  Date:  5  Mar  1995.. ..all  rights  reserved... 


% 


%  Create  the  Efficiency  matrix  (call  efffunc) 

%  initialize  the  all  data  matrices  and  set  default  terrain 
%  to  flat. 


res=10;  %  10  times  size  of  efficiency  matrix 

nmot  =  0;  %  initialize  matrix 

[effmat  rpmin  Tlin  nmot  isnew]=efffunc(res,nmot); 


datamat  =  zeros(3,23); 

%  parameter  initialization 

strmat  =  [' 

hildef=0; 

[distvec  hiUang  initelev  hUlvec  hildef]=hillfunc(hildef); 


%  Initialize  the  user  control  window  and  set  the  position 
% 

fignumber=figure( ... 

’Position',[0  0  570  200],... 

’Color',[001],... 

’NumberTide','off , ... 

’Name','Welcome  to  the  Solar/Electric  Vehicle  Simulator’, ... 
'Resize','off , ... 

’Pointer','arrow'); 


axes('Units','normalized', ... 
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•Position', [.15  .25  .75  .45]); 

plot([0  2  2.2 .2 .4  2.4], [0  0 .1 .1 .2  .2],'m', ... 

[3  3.4  5.4  5.2  3.2],[-.01 .2  .2  .1  .l],'m’, ... 

[8  6  6.2  7.4],[0  0 .1  .l],’m’,[6.2  6.4  8.4], [.1 .2 .2],'m', 
[9.4  9  11],[.21-.01 .22],'m',... 

'LineWidth',12) 
axis([-.5  11.5  -.1 .3]) 
axis  off 


^  ********************  He  ***********  * 

%  Create  the  OPEN  SIMULATION  WINDOW  Button 


spevbut  =  uicontroK ... 
’Style',’pushbutton', ... 

'Units',’pixels', ... 

'Position',[320  10  16030],... 

'String', 'Open  Simulation  Window', ... 
•Visible', 'off, ... 

Enable', 'off, ... 

'Callback','newsim'); 


9fc  %  9fe  sfc  9{c  »|c  9fc  9{e  ^  3{e  *  He  ********************  * 

%  Create  WELCOME  Button 


titll  =  Welcome'; 

welcscm  =['  *  Welcome  to  the  Solar  Vehicle  Simulator. 
'  *  Please  make  your  selections  from  the 
'  *  Control  Buttons  in  the  window  and  follow  ’ 

'  *  the  directions  at  each  prompt  in  the  ' 

'  *  Matlab  Command  window.  ' 

»  He  * 

»  He  * 

»  He  * 

'  *  --  Enjoy  the  ride!  -  ']; 


welcbut  =  uicontrol( ... 

'Style', 'pushbutton', ... 
'Units','pixels', ... 

'Position',[5  10  100  30], ... 
Visible’,'off , ... 

'String', Welcome', ... 

Enable','off, ... 

'Callback','helpsp(titl  1  ,welcscm)'); 


^He*********  ************************* 
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%  Create  the  HELP  Button 
% 

titl2  =  'Solar/Electric  Vehicle  Simulator  Help'; 
helpscm  =  ['  ’ 

'  The  SPEV  Simulator  will  assist  you  in  ' 

'  optimzing  the  configuration  of  your  ' 

'  vehicle  by  providing  user  controlled  ' 

'  parameter  variation  and  continuous  ' 

'  display.  Your  vehicle  characteristics  ' 

'  may  be  easily  changed  without  the  initial' 

'  design  costs.  ' 

I  t 

'  Use  the  "User  Control  Buttons"  to  input  ' 

'  the  vehicle  characteristics,  then  change ' 

'  individual  characteristics  by  following  ' 

'  the  directions  in  the  Matlab  Control  ' 

'  Window. 

I  I 

'  Click  on  "Open  Simulation  Window"  to  open' 

'  the  Simulink  model  for  the  vehicle.  ' 

I  I 

'  Use  the  "Throttle"  to  adjust  the  vehicle ' 

'  speed  and  "Cloud  Cover"  to  adjust  the  ' 

'  current  from  the  solar  array  (simulate  ' 

'  cloud  cover) ...  see  the  results. 

helpbut  =  uicontrolf ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

’Position',[110  10  100  30],... 

'Visible','off , ... 

'String','Help', ... 

'Enable','off , ... 

'Callback','helpsp(titl2,helpscm)'); 

9|c  9|e  s)c  :(c  9|c  %  4: 3(e  :fe  9{e  9ic  *  3fe  ajc sfe  He  3ie  3ie  sK  He  ^  ^ 

%  Create  the  CLOSE  Button ...  closes  the  user  control  window 

closebut  =  uicontrolf ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

'Position',[215  10  100  30], ... 

'Visible','off , ... 

'String', 'Close', ... 

Enable','ofF, ... 
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'CaUbackVclose(gcf)'); 


***********************  ********** 


%  Create  all  User  Control  Buttons 


%  LOAD  INITIAL  DATA ...  loads  all  initial  and  default  data  for 
%  first  run  simulations. 

%  The  default  is  mechanical  data  from  UC  Berkeley’s  solar  car 
%  with  a  Unique  Mobility  1 1.3  Hp  motor,  tracking  array,  and 
%  lead-acid  batteries. 

% 

firstbut  =  uicontrol( ... 

'Style','pushbutton', ... 

'Units','pixels', ... 

’Position',[5  170  140  30], ... 

Visible','off , ... 

'String’,'Load  Initial  Data', ... 

Enable', 'off, ... 

'Callback','[datamat  strmat]=spevfunc("first",datamat,strmat);'); 


%  RUN  TIME ...  controls  the  total  run  time  desired,  the  time 
%  of  day  to  start  simulation,  and  the  integration  time  step  (dt) 

%  which  controls  the  speed  of  simulation. 

% 

runbut  =  uicontrol( ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

’Position’,[150  170  100  30], ... 

’Visible’,’off, ... 

'String',’Run  Time', ... 

Enable', 'off, ... 

'CaUback','[datamatstrmat]=spevfunc("runtime",datamat,strmat);'); 


%  MOTOR  PARAMETERS...  allows  the  user  to  change  the  motor  type  or 
%  individual  parameters  ...  default  is  a  1 1.3Hp/8.4KW  motor. 

% 

motbut  =  uicontrol( ... 

'Style','pushbutton', ... 

'Units','pixels', ... 

'Position',[255  170  145  30], ... 

'Visible','off , ... 

'String','Motor  Parameters', ... 
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'Enable', 'off , ... 

'Callback', '[datamat  strmat]  =spevfunc("motor",datamat,straiat);') ; 


%  VEHICLE  PARAMETERS  ...  allows  the  user  to  change  the  vehicles 
%  mechanical  parameters.  The  default  data  is  obtained  from  the 
%  UC  Berkeley  (CALSOL)  team  vehcile. 

% 

vehbut  =  uicontrol( ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

'Position',[405  170  155  30], ... 

Visible','off , ... 

'String','Vehicle  parameters', ... 

Enable','off , ... 

'Callback', '[datamat  strmat]=spevfunc("vehicle",datamat,strmat);'); 


%  ARRAY  TYPE ...  allows  the  user  to  choose  between  tracking, 
%  flat,  or  curved  array  configurations  and  specify  the  peak  current 
%  expected.  Default  is  tracking  @11  amps  peak 

% 

araybut  =  uicontrol( ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

'Position',[5  135  140  30],... 

Visible',’off , ... 

'String', 'Array  Type', ... 

Enable','off , ... 

'Callback',  '[datamat  strmat]=spevfunc("saray",datamat,strmat);'); 


%  BATTERY  TYPE ...  choose  between  lead-acid  or  silver-zinc  batteries. 
%  battery  voltage  is  set  to  maintain  a  100  volt  bus  voltage. 

%  Default  is  lead-acid  battery. 

% 

battbut  =  uicontroK ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

Eosition',[150  135  100  30], ... 

Visible','off , ... 

'String','Batt  Type', ... 

Enable', 'off , ... 

'Callback' ,'[datamatstrmat]=spevfunc("battery",datamat,strmat);'); 
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%  CHANGE  TERRAIN  PROFILE...  allows  user  the  change  the  driving  terrain. 
%  Choose  between  flat,  single  hill,  or  terrain  profile  from  a  data  file.  Data 
%  must  be  loaded  from  a  2  column  "filename.dat"  text  file.  First  column 
%  is  distance  from  start  in  1/10  mile  increments  and  the  corresponding 
%  elevation  is  feet  above  sea  level.  Default  is  flat  terrain. 

% 

hillbut  =  uicontroK ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

'Position',[255  135  145  30], ... 

Visible','off , ... 

'String','  Change  Terrain  Profile', ... 

'Enable', 'off, ... 

'Callback','[distvec  hdlang  initelev  hiUvec  hildefj=hitlfunc(hildef);'); 


%  SHOW  ALL  DATA ...  shows  all  significant  data  into  simulator. 
%  Useful  after  making  changes. 

% 

showbut  =  uicontroK ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

'Position',[405  135  155  30], ... 

'Visible','off , ... 

'Strmg’,'Show  All  Data', ... 

'Enable', 'off, ... 

'Callback', '[datamat  strmat]=spevfunc("showall",datamat,strmat);'); 


%  MOTOR  EFFICIENCY  PLOT ...  3-dimensional  plot  of  motor  efficiency. 
%  Each  motor  must  have  an  efficiency  matrix. 

% 

effbut  =  uicontroK ... 

'Style',’pushbutton', ... 

'Units','pixels', ... 

'Position',[405  100  155  30], ... 

'Visible','off, ... 

'String',’Motor  Efficiency  Plot', ... 

Enable','off , ... 

'CaIlback','e^ltf(rpmin,Tlin,effmat,isnew);'); 


%  CHANGE  MOTOR  EFFICIENCY ...  change  motor  efficiency. 
%  Each  motor  must  have  an  efficiency  matrix. 
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% 

effmotbut  =  uicontrol( ... 

'Style', 'pushbutton', ... 

'Units','pixels', ... 

'Position',[235  100  165  30], ... 

'Visible','off , ... 

'String','Change  Motor  Efficiency', ... 

'Enable','off , ... 

'Callback','[effmat  rpmin  Hin  nmot  isnew]=efffunc(res,nmot);'); 

%  Turns  on  all  buttons  .... 

% 

set([welcbut  spevbut  helpbut  closebut  firstbut  runbut  motbut ... 
vehbut  araybut  battbut  hillbut  showbut  effbut  effmotbut], ... 
'Enable',’on', ... 

'Visible','on’); 


B.  FUNCTION  SPEVFUNC.M  TO  CHANGE  SYSTEM  PARAMATERS 
%  Solar  Powered  Electric  Vehicle  Simulator. 

%  Developed  by  Steven  J.  Roerig,  LT,  USN,  at  the  Naval-Post  Graduate 
%  School,  Monterey,  Ca.  Date:  7/26/94 
% 

%  This  function  will  change  the  system  paramters  when  called. 

% 


function  [spevdata,strdata]  =  spevfunc(action,data,strmat); 


%  Load  the  initial  data  and  vehicle  parameters... 

% 

if  strcmp(action,'first'); 

dispCSince  this  is  your  initial  run,  the  default  data  has  been  entered') 
dispCfor  a  11.3  Hp,  tracking  array,  Pb-Acid  vehicle'); 
data(l,(l;10))  =  [.5  3600*9  9  24  0 .014 .036  25e-6/.036  20 .184]; 
data(  1,(1 1:23))  =  [.16 .0105  1.2  295.8  .98  35  100  0  0 .26  5  0  0]; 
data(2,:)  =  1.1*[10  10  10  10  10  10  10  10  10  10 ... 

10  10  10  10  10  10  10  9.9  9.8  0  0  0  0]; 
strmat(l,:)  =  'Tracking  '; 

data(3,:)  =  60*([.6 .9  1.25  1.4  1.55  1.65  1.7  1.75  1.8  1.84  1.88  ... 

1.91  1.94  1.96  1.98  2.0  2.02  2.04  2.06  2.07  2.08  0  0]); 
strmat(2,;)  =  'Pb-Acid  ’; 

end 


%  Modify  the  run  time  parameters. 
% 
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if  strcmp(action, 'runtime') 

dispCNote:  All  run  times  will  begin  at  9  AM.  Maximum  run  time') 
dispCis  8  Hours.') 

tmin  =  inputC  Enter  length  of  simulation  run  (minutes): '); 
data(l,2)  =  60*tmin;  %final  time  in  seconds 

timein  =  inputC  Simulation  start  time  (btwn  0900  and  1800): '); 
data(l,3)  =  timein/100; 

data(l,l)  =  inputC  Time  step  (how  slow/fast  to  run:  btwn  .0001  and  1  sec): '); 
end 


%  Modify  the  motor  characteristics 
% 

if  strcmp(action,'motor') 
dispCMotor  Characteristics  -') 
data(l,4)  =  input('  Number  of  poles  (P): '); 
data(l,5)  =  inputC  Damping  coefficient  (B): '); 
data(l,6)  =  input('  Inertia  (J): '); 
data(l,7)  =  inputC  Winding  resistance  (ra): '); 
data(l,8)  =  inputC  Electrical  time-constant  (ta): '); 
data(l,9)  =  inputC  Back-EMF  constant  (ke): '); 
data(l,10)  =  inputC  Torque  constant  (kt): '); 
end 


%  Modify  the  vehicle  characteristics 
% 

if  strcmp(action,’vehicle') 
dispCVehicle  Characteristics  -') 
data(l,l  1)  =  inputC  Drag  coefficient  (Cd): '); 
data(l,12)  =  inputC  Rolling  resistance  coefficient  (Cr):  *); 
data(  1,13)  =  inputC  Frontal  area  (A  in  m'^2): '); 
data(l,  14)  =  inputC  Vehicle  mass  (Mv  in  kg): '); 
data(l,20)  =  inputC  Drive  wheel  radius  (r  in  meters): '); 
data(l,21)  =  inputC  Drive  ratio  (belt,  chain,  or  gear  -  1:N)  (N=): '); 
effm  =  inputC  Drive  train  efficiency  (%): '); 
data(l,15)  =  effm/100; 
end 


%  Change  solar  array  type  and  peak  output  current 
% 

arr.dat  =  [10  10  10  10  10  10  10  10  10  10 ... 

10  10  10  10  10  10  10  9.9  9.8  0  0  0  0]; 
if  strcmp(action,'saray') 
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ar_choiee  =  input('Solar  Array  type:  Flat(l),Curved(2),Tracking(3):  Number  - '); 
peak  =  input('What  is  your  expected  peak  current: '); 
if  ar_choice  =  1 

data(2,:)  =  (peak/10)* [.72  .78  .83  .88  .92  .97  .98  .98  .99  .985  .98  ... 

.95  .92  .90 .85  .82 .77  .67  .48  0  0  0  0].*arr_dat; 
strmat(l,:)  =  'Flat 
elseif  ar_choice  ==  2 

data(2,:)  =  (peak/10)*[.59  .63  .68  .75  .79  .85  .87  .87  .88  .87  ... 

.86  .85  .82  .78  .72  .69  .64 .55  .48  0  0  0  0].*arr_dat; 
strmat(l,:)  =  'Curved 
elseif  ar_choice  ==  3 
data(2,:)  =  (peak/10)*arr_dat; 
strmat(l,:)  =  'Tracking  '; 

end 


end 


%  Change  battery  type. 

% 

if  strcmp(action, 'battery') 

bat_choice  =  input('Battery  type:  Pb-Acid(l),Ag-Zn(2):  Number  - '); 
if  bat_choice  ==  1 

data(3,:)  =  60*([.6 .9  1.25  1.4  1.55  1.65  1.7  1.75  1.8  1.84  1.88  ... 
1.91  1.94  1.96  1.98  2.0  2.02  2.04  2.06  2.07  2.08  0  0]); 
strmat(2,:)  =  'Pb-Acid  '; 

elseif  bat_choice  ==  2 

data(3,:)  =  68*([.8  1.35  1.45  1.47  1.48  1.49  1.495  1.497  ... 

1.5  1.5  1.5  1.5  1.5  1.5  1.5  1.5  1.55  1.7  1.8  1.85  1.9  0  0]); 
strmat(2,:)  =  'Ag-Zn 
end 

data(l,16)  =  inputCWhat  is  the  A-h  rating  of  your  battery: '); 
data(l,17)  =  inputCWhat  is  your  initial  state  of  charge  (percent): '); 
end 


%  Show  all  data  upon  initialization  and  whenever  called. 

% 

if  strcmp(action,'showall')  I  strcmp(action, 'first') 
dispC) 

dispCThe  below  information  is  currently  is  loaded  into  the  simulator:') 
disp(['  Run  time:  ',num2str(data(l,2)/60),'  minutes']) 
disp(['Start  time:  ',num2str(data(l,3))]) 
disp(['  Time  step:  ',num2str(data(l,l))]) 


dispC) 
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dispCMotor  Characteristics ') 
disp(['  Number  of  poles:  ',num2str(data(l,4))]); 
disp(['  Damping  coefficient:  ',num2str(data(l,5))]); 
disp(['  Inertia:  ',num2str(data(l,6))]); 
disp(['  Armature  resistance:  ',num2str(data(l,7))]); 
disp(['  Electrical  time-constant:  ',num2str(data(l,8))]); 
disp(['  Back-EMF  constant:  ',num2str(data(l,9))]); 
disp(['  Torque  constant:  ',num2str(data(l,10))]); 

dispO; 

dispCVehicle  Charcteristics '); 

disp(['  Drag  coefficient:  ',num2str(data(l,l  1))]); 

disp(['  Rolling  resistance  coefficient:  ’,num2str(data(l,12))]); 

disp(['  Frontal  area  (meters  sqrd):  ',num2str(data(l,13))]); 

disp(['  Vehicle  mass  (kilo-gram):  ',num2str(data(l,14))]); 

disp(['  Drive  wheel  radius  (meters):  ',num2str(data(l,20))]) 

disp([’  Drive  ratio  (N):  l:',num2str(data(l,21))]); 

disp(['  Drive  train  efficiency:  ',num2str(data(l,15)*100), '%']); 

dispC’); 

disp(['Solar  Array  Type:  ’,strmat(l,l:21)]); 

dispO; 

disp(['Battery  Type:  ',strmat(2,l:21)]); 


dispC) 


dispCIf  you  need  to  make  changes  in  your  data,  just  click  on  the') 
dispCitem  in  the  Menu  Bar.') 
end 


%  Output  to  the  system. 

% 

spevdata  =  [data(l,l:23);data(2,l:23);data(3,l:23)]; 
strdata  =  [stnnat(l,l:23);strmat(2,l:23)]; 

%  end  of  function 


C.  FUNCnON  HILLFUNC.M  TO  CHANGE  TERRAIN  PROFILE 

%  This  function  performs  all  hill  and  terrain  profile  modifications 
%  when  called  from  the  user  control  window. 

% 


function  [distvecjiillang,initelev,hillvec,hildef]=hillfunc(hildef); 
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ifhildef~=0 

whichone  =  inputCFor  a  single  hill  enter  - 1,  hill  profile  from  file  enter  -  2,  or  flat  driving 
enter  -  3: '); 

dispC) 

else 

whichone  =  3; 
end 


%  For  a  single  hill  terrain  profile. 

% 

if  whichone  ==  1 

stl  =  input(’How  many  miles  is  the  hill  from  simulation  start? '); 

st2  =  inputCHow  long  is  the  hill  (in  miles)? '); 

ell  =  inputCWhat  is  your  initial  evelation  (feet)? '); 

el2  =  inputCWhat  is  will  be  your  final  elevation? '); 

distvec  =  0:.l;stl+st2; 

if(el2-ell)>0 

hillvec  =  [ell*ones(size(l:length(0:.l:(stl-.l))))  ell:((el2-ell)/(10*st2)):el2]; 
elseif  (el2-ell)  <0 

hiUvec  =  [ell*ones(size(l:length(0;.l:(stl-.l))))  ell:-((el2-ell)/(10*st2)):el2]; 
end 

type  =  'Single  HiU'; 

%  or  for  a  flat  driving  profile 
% 

elseif  whichone  ==  3 
ifhildef~=0 

stl  =  inputCHow  many  miles  will  you  be  driving? '); 
ell  =  inputCAt  what  elevation  will  you  be  driving  (feet)? "); 
else 

stl  =  20; 
ell  =  100; 
end 

distvec  =  0:.l:stl; 
ifell~=0 

hillvec  =  [ell*ones(size(l:length(0:.l:stl)))]; 
else 

hillvec  =  [2eros(size(l:length(0:.l:stl)))]; 
end 

type  =  Tlat  Terrain'; 
hildef=l; 

%  or  if  a  terrain  profile  is  to  be  loaded  from  a  text  file. 

% 
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elseif  whichone  =  2 

dispCNOTE;  The.  input  file  must  be  set  up  with  two  column  vectors.  The  first'); 
dispCcolumn  contains  the  distance  from  start  of  simulation  in  tenths  of  miles,'); 
dispCand  the  second  column  contains  the  elevation  above  sea  level  in  feet.  The'); 
dispCfile  must  be  a  text  fiile  of  the  form  "filename.dat".'); 
dispC) 

filename  =  input('Input  the  filename:  ','s'); 
fid=fopen(filename); 
dat  =  fscanf(fid,'%g  %g',[2inf]); 
fclose(fid); 
distvec=  dat(l,:); 
hillvec=  dat(2,:); 
type  =  'Hilly  Terrain'; 
end 

%  compute  the  hill  angle  (alpha)  and  output  to  system 

% 

initelev=hillvec(  1 , 1 ) ; 
for  i=2:length(hillvec) 

hillang(i-  l)=asin((hillvec(i)-hillvec(i- 1))/528); 
end 

hiUang(length(hitlvec))=hillang(length(hillang)); 

%  Plot  the  updated  terrain  profile  for  use  during  simulation 
% 

Figures  =  get(0, 'Children'); 
for  i  =  l:length(Figures) 

if  strcmp(get(Figures(i),'Type'),  'figure') 
if  strcmp(get(Figures(i),  'Name'),  Terrain  Profile') 
close(Figures(i)) 
break; 
end 
end 
end 

figure('NumberTitle','ofF,'Name',Terrain  Profile', ... 

'NextPlot','add’,'Position',[0  0  450  300]); 
plot(distvec,hitlvec),grid,ylabel('feet  above  sea  level') 

axis([0  max(distvec)  (min(hillvec)-(.l*min(hillvec)))  (max(hillvec)+(.l*max(hitlvec)))]) 
xlabel('distance  in  miles'),title(type) 

D.  FUNCTION  EFFFUNC.M  TO  EXPAND  SIZE  OF  EFFICIENCY  MATRIX 
%  This  function  will  expand  the  efficiency  matrix  in  order 
%  to  decrease  interpolation  errors. 

% 
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function  [efficiency ,rpmin,Tlin,nmot,isnew]  =  efffunc(N,nmot); 

%  Efficiency  for  the  1 1 .3  Hp  motor 
% 


isnew=  127; 
if  nmot  ==  1 

isnew=input('Is  this  the  127  or  86  Unique  motor?  ’); 
end 


if  isnew  =  86 

effmat=[35  55  67  68  69  70  71  70  63  55; 

58  68  72  74  75  74  72  70  68  66; 

68  74  78  78.5  78.6  76  74  72  70  69; 

72  77  82  82  81.9  81  80  77  73  70; 

74  79  84  84.4  85  84.8  83  82  80  78; 

78  80  86  87.3  87.2  87.1  86.5  84.5  83  81.5; 

81  84  89  89.7  90.1  89.8  88  87  86  85.5; 

82  86  90.5  92  92.2  92.8  91  89.2  88.5  87; 

83  87  91.7  94  93.5  92.5  91.6  89.8  88.6  87.1; 

84  87.5  93  94.5  94.2  93  91.9  90  89.5  88.5]; 

elseif  isnew  ==  127 

effmat=[10  20  32  39  38  40  51  41  28  18; 

18  36  60  71  70  69  72  63  55  45; 

20  50  76  85.5  82  80.5  79  75  71  64; 

30  60  85.5  86  85.5  84.8  83  82  79  74; 

33  71  86  87.5  86.5  86.5  85.5  85  83.5  81; 

38  80  87  88.5  89  88  87  86  95.5  84; 

40  84  89  90  90.5  90  89  88  87  86; 

50  85  89.5  91  92  91.5  91  90  90.5  89.5; 
50  86  90  92  93  93.5  91  91.5  90.5  89.5; 
48  86.5  90.5  92.5  94  93.5  93  92  91  90.5]; 
nmot=  1; 
end 


[n  m]=size(ef&nat); 


%  Increases  the  size  of  the  matrix  from  n  x  m  to  N*n  x  N*m 
%  with  linspace  interpolation. 

% 

forj=l:n 
for  i=2:m 

eftemp(i-l,:)=linspace(effmat(j,i-l),effmat(j,i),N); 
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end 

[c  d]  =  size(eftemp); 
for  a=l  :c 
for  b=l:d-l 

neweffmat(j,(((N-l)*(a-l))+b))  =  eftemp(a,b); 
end 
end 
end 


new2=neweffmat'; 
[n  m]=size(new2); 


forj=l:n 
for  i=2:m 

eftemp2(i- 1  ,:)=linspace(new2(j,i-  l),new2(j,i),N)5 
end 

[c  d]  =  size(eftemp2); 
for  a=l  :c 
forb=l:d-l 

neweffmat2(j,(((N-l)*(a-l))+b))  =  eftemp2(a,b); 
end 
end 
end 

efficiency  =  neweffmat2'; 

[n  m]=size(efficiency); 

%  Corresponding  torque  and  rpm  vectors. 

% 

ifisnew==127  %  Unique  DR127/CR10 
rpmin=linspace(25 ,5(X)0,n); 

Tlin=linspace(.  l,20.25,m); 
elseif  isnew  =  86  %  Unique  DR086/CR10 
rpmin=linspace(200,5600,n); 

Tlin=linspace(.4, 15.2, m); 
end 


E.  FUNCTION  EFFPLTF.M  FOR  SURFACE  PLOT  OF  EFFICIENCY  MATRIX 

%  This  function  performs  a  surface  plot  of  the  motor 
%  efficiency 

function  effpltf(rpmin,Tlin,effmat,isnew); 
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figure('NumberTide','ofF, 'Name', 'Efficiency  Window') 
colormapOet) 

surf(Tlin,rpmin,effmat),grid 

axis([0  max(Tlin)  0  max(rpmin)  0  100]) 

shading  flat 

colorbar 

xlabel('Torque'), 

ylabel('RPM'), 

zlabel('Eff(%)'), 

if  isnew  =127 

title('Motor  Efficiency  for  Unique  DR127/CR10'); 
elseif  isnew  ==  86 

title('Motor  Efficiency  for  Unique  DR086/CR10’); 
end 


F.  FUNCTION  SPEVAN2.M  FOR  ON-LINE  DISPLAY  WINDOW 

function  [sys,x0]=spevan2(t,x,u,flag,fignamel); 

%  spevan2  S-fUnction  for  animating  the  on-line  display  of  dynamic 
%  output.  To  see  the  specific  input  corresponding  to  the  parameter 
%  u[],  open  the  "On-Line  Display"  block  in  the  simulation  window. 

%  Steven  Roerig.  7/25/9 

global  HNDLl 
if  abs(flag)  =  2 

if  (x(l)  ~=  Inf)  %  screen  has  been  initialized 
hndls=get(x(13),'UserData’); 

if  u(5)  >  -.01  &  u(5)  <  .01  %  set  output  to  zero 

set(hndls(  1), 'String', num2str(0.00)); 
else 

set(hndls(  1 ) , 'String', num2str(u(5))) ; 
end 


set(hndls(2),’String',num2str(u(3))); 

set(hndls(3),'String',num2str(u(4))); 

if  u(l)  >  -.001  &  u(l)  <  .001  %  set  output  to  zero 

set(hndls(4),'String',num2str(0.00)); 
else 

set(hndls(4),'String',num2str(u(  1))); 
end 

set(hndls(5),'String',num2str(u(2))); 
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set(hndls(8),'Strmg',num2str(u(7))); 
set(hndls(12),’String',num2str(u(6))); 
set(hndls(9), 'String', num2str(u(8))); 

if  u(9)  >  -.001  &  u(9)  <  .001  %  set  vehicle  "stopped" 
set(hndls(  1 0),'S  tring',num2str(0.00)) ; 
set(hndls(ll), 'String', 'Stopped.', 'Color', [1  1  0]); 
elseif  u(9)  <  -.001  %  set  vehcile  "backward" 

set(hndls(10),'String',num2str(u(9))); 
set(hndls(ll),'String','Going  Backward!  !','Color',[l  0  0]) 
else  %  set  vehicle  "going  forward" 

set(hndls(10),'String',num2str(u(9))); 
set(hndls(ll),'String','GoingForward.','Color',[0  1  1]) 

end 

if  u(5)  >  .001  %  set  drawing  current 

set(hndls(6),'XData',l,'YData',2.5,'Color',[0  1  0]); 
set(hndls(7),'String','using'); 
elseif  u(5)  <  -.001  %  set  regenerating  current 

set(hndls(6),'XData',l,'YData',2.5,'Color',[l  0  0]); 
set(hndls(7),'String','regen'); 
else  %  set  motor  stopped 

set(hndls(6),'XData',l,'YData',2.5,'Color',[0  0  0]); 
set(hndls(7),’String',' '); 
end 

set(hndls(13),'String',num2str(u(10))); 

if  u(l  1)  <=  .01  %  set  output  to  zero 

set(hndls(14),'String',num2str(0.00)); 
else 

set(hndls(14),'String',num2str(u(ll))); 

end 

%  output  the  speed  to  analog  display 
y2=[1.5  1.5+1.  l*sin((2*pi*u(9)/60))]; 
x2=[10.5  10.5-l.l*cos((2*pi*u(9)/60))]; 
set(hndls(15),'XData',x2,'YData',y2); 
set(hndls(16),'XData',10.5,'YData',1.5); 
set(hndls(17),'String',num2str(u(12))); 
drawnow; 
end 

sys=[u;x(  13)];  %  initialize  system  parameters 

elseif  flag  =  0  %  Initialize  the  figure  "display"  window 
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sys=[0;12+l;0;12;0;0]; 

Figures  =  get(0,'Chil'); 
yesfig  =  0; 

for  i  =  l:length(Figures) 
if  strcmp(get(Figures(i),  'Name'),  'Output') 
yesfig  =  1; 
whichone  =  i; 
end 
end 

if  yesfig  ==  1 

curhndl  =  Figures(whichone); 
x0=[lnf;0;0;0;0;0;0;0;0;0;0;0;curhndl]; 
set(curhndl,'Userdata',HNDLl); 
else 

curhndl  =  animl(fignamel);  %  draw  the  speed  dial 
cir  =  pi:-.035:-pi; 
for  i=l:length(cir) 
xl(i)  =  10.5  +  1.2*cos(cir(i)); 
yl(i)  =  1,5  +  1.2*sin(cir(i)); 
end 

for  i=  1  :length(cir)/ 1 5 
xla(i)  =  xl(i*length(cir)/12); 
yla(i)  =  yl(i*length(cir)/12); 
end 

x2=[10.5  9.4]; 
y2=[1.5  1.5]; 

hold  on  %  everything  below  initializes  the  display  window 
plot(x  1  ,y  1 ,  'r')  ,plot(x  1  a,y  1  a,',  w',  'MarkerSize',  11); 
axis([0120  3.8]); 
axis  off 

text(8.9,1.55,'0','FontName’,Times  New  Roman','FontSize',[10]) 
text(10.3,3.,'15','FontName', 'Times  New  Roman’,'FontSize',[10]) 
text(l  1.9, 1.55, '30', 'FontName', 'Times  New  Roman','FontSize',[10]) 
text(10.4,0,'45','FontName','Times  New  Roman’,'FontSize’,[10]) 
text(10.1,3.5,'Mph’,'FontName',Times  New  Roman',’FontSize',[10]) 
text(3.35,3,’Array  Current  (amps):','FontName',Times  New  Roman','FontSize',[10]) 
text(3.95,2.5,'Run  time  (secs):','FontName',Times  New  Roman','FontSize’,[10]) 
text(3.72,2,'Odometer  (miles):', 'FontName',Times  New  Roman','FontSize',[10]) 
text(3.5,1.5,'Speedometer  (mph);','FontName','Times  New  Roman', 'FontSize',[10]); 
text(-.l,3,'  Motor  Current  (amps):','FontName',Times  New  Roman','FontSize',[10]) 
text(.92,2,'Ahrs  Used:','FontName’, 'Times  NewRoman','FontSize',[10]) 
text(l,1.5,'SOC  (%):', 'FontName’,'Times  New  Roman’,'FontSize',[10]) 
text(.05,.5,'  Vehicle  Load  (N-m):’,'FontName','Times  New  Roman', 'FontSize',[10]); 
text(7.3,3,'  RPM;','FontName', 'Times  New  Roman','FontSize',[10]); 
text(6.4,2.5,'Hill  Grade  (%):','FontName','Times  New  Roman','FontSize',[10]); 


text(-.  1,1, 'Battery  Voltage  (volts):', 'FontName','Times  New  Roman’,'FontSize',[10]); 
text(3.35,.5,'Motor  Efficiency  (%):','FontName','Times  New  Roman','FontSize',[10]); 
HNDLl(l)=text(2.2,3,'  ',’EraseMode','background','FontName','Times  New 
Roman’,'FontSize',[  10]); 

HNDLl(2)=text(2.2,2,'  ','EraseMode','background','FontName',’Times  New 
Roman','FontSize',[  10]); 

HNDL 1  (3)=text(2.2, 1.5,'  ','EraseMode','background','FontName','Times  New 
Roman','FontSize',[  10]); 

HNDL  1  (4)=text(2.2,.5,'  ','EraseMode','background','FontName', 'Times  New 
Roman','FontSize',[  10]); 

HNDLl(5)=text(5.6,3,'  ','EraseMode’,'background','FontName','Times  New 
Roman','FontSize',[  10]); 

HNDLl(6)=plot(l,2.5,'.','MarkerSize',30,'EraseMode','background'); 

HNDL  1  (7)=text(  1 .2,2.5,'  ’,'EraseMode','background','FontName','Times  New 
Roman', 'FontSize',[  10]); 

HNDLl(8)=text(5.6,2.5,'  ',’EraseMode','background','FontName','Times  New 
Roman','FontSize',[  10]); 

HNDLl(9)=text(5.6,2,'  ','EraseMode','background','FontName','Times  New 
Roman','FontSize',[  10]); 

HNDLl(10)=text(5.6,1.5,'  ','EraseMode','background','FontName','Times  New 
Roman’,'FontSize’,[  10]); 

HNDLl(ll)=text(4,l,'  ','Color',[l  0  0],EraseMode','background','FontName','Times 
New  Roman','FontSize',[10]); 

HNDLl(12)=text(2.2,l,'  ’,'EraseMode','background','FontName','Times  New 
Roman', 'FontSize',[  10]); 

HNDLl(13)=text(5.6,.5,'  ','EraseMode','background',’FontName', Times  New 
Roman', ’FontSize',[  10]); 

HNDLl(14)=text(8,3,'  ','EraseMode','background','FontName','Times  New 
Roman','FontSize',[  10]); 

HNDLl(15)=plot(x2,y2,'c','LineWidth',1.5,'EraseMode','background'); 
HNDLl(16)=plot(10.5,1.5,'.','MarkerSize',20,'EraseMode','background'); 
HNDLl(17)=text(8,2.5,'  ’,'EraseMode','background','FontName',Times  New 
Roman', 'FontSize',[  10]); 
hold  off 

x0=[lnf;0;0;0;0;0;0;0;0;0;0;0;curhndl]; 

set(curhndl,'UserData',HNDLl); 

end 

else 

sys=n; 

end; 


G.  FUNCTION  ANIM1.M  TO  INITIALIZE  ON-LINE  DISPLAY  WINDOW 

function  [curhndl]  =  animl(fignamel); 

%animl  Initializes  a  figure  for  current  animations. 
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%  Steven  Roerig.  25  July  1994,  for  SPEV 


%  Now  initialize  the  whole  figure... 
position=[5  5  825  250]; 
curhndl=figure( ... 
'Name',fignamel, ... 
'NumberTitle’.'ofF, ... 
’BackingStore','off, ... 
Position’, position); 
axes( ... 

'Units','normalized', ... 
'Position',[.05  .1 .9  .75], ... 
'Visible’,'on', ... 
DrawMode','fast'); 


cla  reset 
axis  off 


H.  FUNCTION  PLOTS.M  TO  GENERATE  FINAL  OUTPUT  PLOTS 

%  This  routine  generates  the  plots  for  specfice  outputs  as 
%  well  as  averages  for  the  entire  simulation. 

%  Plot  the  miles  per  hour  and  the  total  force  on  the  vehicle 
% 

figure(Position’,[0  340  300  300]) 
subplot(2, 1 ,1  ),plot(miles,mph),grid 
title('miles  per  hour') 

subplot(2,  l,2),plot(miles,Ftot),grid,xlabel('miles  driven') 
title('Total  force  on  the  vehicle  (N)') 

%  Plot  the  current  drawn  by  the  motor  and  the  total  amp  hours  used 

% 

figure('Position',[0  0  300  300]) 
subplot(2,l,l),plot(mUes,iq),grid, 
title('current  drawn  by  motor') 

subplot(2, 1 ,2),plot(miles,Ahours),grid,xlabel('miles  driven') 
title('Amp-hours  used') 

%  Plot  the  power  delivered  to  the  drive  wheel  from  the  motor  and 
%  the  power  required  at  the  rear  wheel  as  a  function  of  the  forces 
% 

figure('Position',[340  0  300  300]) 
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subplot(2, 1 ,  l),plot(miles,pwrdel),grid, 
title('Power  delivered  to  the  rear  wheels') 

subplot(2,l,2),plot(miles,pwrreq./((datamat(l,15)/100)*effout)),grid,xlabel('miles  driven') 
title('Power  required  at  the  rear  wheels') 

%  Plot  the  efficiency  of  the  motor  for  the  simulation 
% 

figum('Position',[340  340  400  300]) 
subplot(2, 1 , 1  ),plot(miles  ,effout),grid, 
title('Motor  efficiency') 

%  Display  the  averages  beyond  the  transient  system  startup 
% 

subplot(2,l,2) 
axis([0  2  0  2]), axis  off 
text(.22,1.5,'Average  Speed: 

','FontSize',[9]),text(.8,1.5,num2str(mean(mph(100:length(mph)))),'FontSize',[9]) 
text(  1 , 1 .5  ,'mph','FontSize',[9]) 
text(0,l, 'Average  Current  Drawn: 

','FontSize',[9]),text(.8,l,num2str(mean(iq(100:length(iq)))),'FontSize',[9]) 
text(  1 , 1  ,'Amps','FontSize',[9]) 
text(-.05,.5,'Average  Motor  Efficiency: 

’,'FontSize',[9]),text(.8,.5,num2str(mean(effout(100:length(effout)))),'FontSize’,[9]) 

text(l,.5,'%','FontSize',[9]) 

text(1.3,1.5,'Average  Load: 

’,'FontSize',[9]),text(l  .8, 1 .5aium2str(mean(n2(100:length(n2))));FontSize',[9]) 

text(2,1.5,'Nm’,'FontSize',[9]) 

text(1.34,l,'Miles  Driven: 

’,'FontSize',[9]),text(1.8,l,num2str(miles(length(miles))),'FontSize',[9]) 
text(2, 1  ,'miles','FontSize',  [9]) 
text(1.25, .5, 'Total  Ahrs  used: 

','FontSize',[9]),text(1.8,.5,num2str(mean(Ahours(100:length(Ahours)))),'FontSize',[9]) 

text(2,.5,'Ahrs',’FontSize',[9]) 
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